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ABSTRACT 

Magnetic field-lines in astrophysical plasmas are expected to be frozen-in at scales larger than the ion 
gyroradius. The rapid reconnection of magnetic flux structures with dimensions vastly larger than the 
gyroradius requires a breakdown in the standard Alfven flux-freezing law. We attribute this breakdown 
to ubiquitous MHD plasma turbulence with power-law scaling ranges of velocity and magnetic energy 
spectra. Lagrangian particle trajectories in such environments become "spontaneously stochastic", so 
that infinitely-many magnetic field-lines are advected to each point and must be averaged to obtain 
the resultant magnetic field. The relative distance between initial magnetic field lines which arrive 
to the same final point depends upon the properties of two-particle turbulent dispersion. We develop 
predictions based on the phenomenological Goldreich & Sridhar theory of strong MHD turbulence and 
on weak MHD turbulence theory. We recover the predictions of the Lazarian & Vishniac theory for the 
reconnection rate of large-scale magnetic structures. Lazarian & Vishniac also invoked "spontaneous 
stochasticity" , but of the field-lines rather than of the Lagrangian trajectories. More recent theories 
of fast magnetic reconnection appeal to microscopic plasma processes that lead to additional terms 
in the generalized Ohm's law, such as the collisionless Hall term. We estimate quantitatively the 
effect of such processes on the inertial-range turbulence dynamics and find them to be negligible in 
most astrophysical environments. For example, the predictions of the Lazarian- Vishniac theory are 
unchanged in Hall MHD turbulence with an extended inert ial range, whenever the ion skin depth 5i 
is much smaller than the turbulent integral length or injection-scale Li. 
Subject headings: galaxies: magnetic fields - methods: theoretical - MHD - turbulence 



1. INTRODUCTION 

It is generally believed that magnetic field embedded 
in a highly conductive fiuid preserves its topolo gy for 
all tinie due to ma gnetic fields being frozen-in ( [Alfven 
1942 Parker 1979). Although ionized astrophysical ob- 
jects, like stars and galactic disks, are almost perfectly 
conducting, they show indications of changes in topol- 
ogy, "magnetic reconnection", on dynamical time scales 



( |Parker|197Q] |Lovelace|[T976l [Priest fc Forbes|2QQ2| . Re- 

connection can be observed directly in the solar corona 



(jlnnes et al."1997 Yokoyama & Shibata 



et al. 



tence o: 



T994), 



1995 



Masuda 
the exis 



but can also be interred 
large-sca l e dynamo activity i nside stellar interi 
ors ([Parker''l993| lOssendrijver 2003). Solar fiares usu- 
ally (Sturrock 1966) and 7-ray bursts often ([Z han^j^ 



Galama et al. |1998p are as 



Yan||2011t jFox et al.n2005 _ 
sociated with magnetic reconnection. Mucl 

|eyink@jhu.edu| 



recent work 



has concentrated on showing how reconnection c an be 
rapid in plasmas with very small collisional rates ([Sha^ 



e^al. 1998; Drake 2001; Drake et al. 2006; Daughton 
e^ah 2006), which substantially constrains astrophysics 



applications of the corresponding reconnection models. 

Reconnection occurs rapidly in computer simulations 
due to the high values of resistivity (or numerical re- 
sistivity) that are employed at the low resolutions cur- 
rently achievable. Therefore, if there are situations where 
magnetic fields reconnect slowly, numerical simulations 
do not adequately reproduce astrophysical reality. This 
means that if collisionless reconnection is the only way to 
make reconnection rapid, then the numerical simulations 
of many astrophysical processes, including those in inter- 
stellar media (ISM), which is collisional, are in error. At 
the same time, it is not possible just to claim that the 
reconnection must always be rapid empirically, as solar 
fiares require periods of fiux accumulation time, which 
correspond to slow reconnection. 
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To understand the difference between reconnection in 
astrophysical situations and in numerical simulations, 
one should recall that the dimensionless combination 
that controls the reconnection rate is the Lundquist num- 
be]{^ defined as 5 = L^va/^^ where is the length of 
the reconnection layer, va is the Alfven velocity, and 
A = r]c^ /4:7T is Ohmic diffusivity. Because of the huge as- 
trophysical length-scales involved, the astrophysical 
Lundquist numbers are also huge, e.g. for the ISM they 
are about 10^^, while present-day MHD simulations cor- 
respond to 5* < 10^. As the simulation costs scale as L^, 
where is the size of the box, it is feasible neither at 
present nor in the foreseeable future to have simulations 
with sufficiently high Lundquist numbers. 

Plasma flows at such high Lundquist numbers are 
generically turbulent, since laminar flows are then prey 
to numerous linear and finite- amplitude instabilities. In- 
deed, turbulence is ubiquitous in astrophysical plasmas. 
This is sometimes driven turbulence due to an external 



energy source, such as supern ova in the ISM (Norman 
fc Ferrara||1996[ [FerrierepQQl ), merger events and AG^ 



outflows Tn the intercluster medi um (ICM) ( Subrama- 
nian et al.|j2QQ6| [EnBlin & Vogt 2QQ6| [Chandra n 2005), 



and baroclmic forcing behind shock waves m mterstellar 
clouds. In other cases, the turbulence is spontaneous, 
with available energy released by a ric h array of instabil- 
ities, such as MRI in accretion disks ( Balbus fc Hawleyj 
|T998), kink instab ility of twisted flux tubes in the solar 
corona (Galsgaard fc Nordlund 1997a; Gerrard fc HooclJ 
[2003 ) , etc. Whatever its origin, the signatures of plasma 
turbulence are seen throughout the universe. Turbulent 
cascade of energy leads to long "inertial ranges" with 
power-law s pectra that are widely observed, e.g. in the 
solar win d p^amon et al.| |1998; Bal e et aL]|2005[), in 



1995 ; Chepurnov fc Lazar-J 

and in the IGM "(IS^ecker et al. 2004; Vogt 



the ISM (Armstrong et al. 
Tan 2010) 



fc Enmin |2QQ5| ). Often inertial range spectra cannot be 
directly observed, but only large-scale turbulent fluctu- 
ations, e.g. through their Doppler broadening of line 
spectra. Nevertheless, the power-law ranges are univer- 
sal features of high- Reynolds- number turbulence which, 
even when not seen, can be inferred to be present from 



enhanced rates of dissipation and mixing (|Eyink||2008D. 
Any theory of astrophysical reconnection rnust take into 
account the pre-existing turbulent environment. Even if 
the plasma flow is initially laminar, kinetic energy release 
by reconnection due to some slower plasma process may 
generate vigorous turbulent motion. 

But can turbulen c e itse lf accelerate reconnection? 
Lazarian fc Vishniac ( 1999| ) [henceforth LV99] proposed 
a model of fast reconnection in the presence of sub- 
Alfvenic turbulence in magnetized plasmas. They iden- 
tified stochastic wandering of magnetic field-lines as the 
most critical property of MHD turbulence which permits 
fast reconnection. As we discuss more below, this line- 
wandering widens the outflow region and alleviates the 
controlling constraint of mass conservation. The LV99 
model has been successfully tested recently in Kowal et 
al. (2009^ (see also higher resolution results in ,Lazarian| 



^ The magnetic Reynolds number, whi ch is the r atio of the mag-' 
netic field decay time to the eddy tur nover time, is defined using 
tl ie injection velocity vi as a characterTstic speed instead of the | 



et al. ( 2010| )). The model is radically different from its 
predecessors which also appealed to the effects of turbu- 
lence (see more comp ariso ns in section 5.4). For insta nce, 
unlike Speiserj ( |197QD and [Jacobson fc Moses| ( |1984D the 
model does not appeal to changes of microscopic proper- 
ties of plas ma. The nearest progenitor to LV99 was the 
work of Matthaeus fc Lamkin ( 1985 , 1986 ), who studied 
the problem numerically in 2D MHD and who suggested 
that magnetic reconnection may be fast due to a number 
of turbulence effects, e.g. multiple X point s and turbu- 
lent EMF. However, Matth aeus fc La mkin | (|1985l |1986D 
did not observe the important role of magnetic field-line 
wandering, and did not obtain a quantitative prediction 
for the reconnection rate, as did LV99. 

The success of LV99 in identifying an MHD turbulence 
mechanism for fast reconnection leads, however, to a con- 
flict with certain conventional beliefs. As the predicted 
reconnection velocity is independent of magnetic diffu- 
sivity A, the LV99 theory implies that field-line topology 
should change in MHD plasmas at a finite rate even in 
the limit of infinite Lundquist number. This contradicts 
the accepted wisdom that magnetic field-lines should be 
"nearly" frozen-in to very high-conductivity MHD plas- 
mas. It is implicit in the LV99 theory that the standard 
Alfven Theorem on magnetic-flux conservation must be 
violated for A ^ (Vishniac fc Lazarian |[l999l). 



In a separate, more recent development, Eyink 
has critically examined the concept of froz en-in tie. 



(2011) 
(times 



for tu r bule nt MHD plasmas. See also Eyink fc Aluie 
( |2006D ; [Efink , ( ,2007 2009 ) . A first contribution of that 
work was to establish exact 



Alfven speed va, which is taken in the Lundquist number. 



flux-conservation properties 
for resistive MHD solutions. The magnetic field at a 
point was shown to result from averaging over an infi- 
nite ensemble of line-vectors stochastically advected to 
that point by the fluid velocity perturbed with a Brow- 
nian motion proportional to a/A- For smooth, laminar 
flows the contribution of the random Brownian term 
disappears as A ^ and standard flux-freezing is re- 
covered. However, it was shown in (Eyink | (2011) that 
flux-conservation may remain stochastic even at infinite 
conductivity for rough, turbulent velocities. This re- 
markable behavior is due to "spontaneous stochasticity" 
of Lagrangian fluid particles undergoin g turbulent two- 
article or Richardson diffusion: see Bern ard et al. 
1998D: [Gawedzki fc Vergassola^ (2000); E fc vand en Er 
jnHen (2000, 2001); Chave s eTal. I (|2003), and, ToFa 
review, Ku piainen (20 03). 'I'he wort of Eyink removes 
the objection to LVyS' based on the Alfven theorem for 
high conductivity plasmas, since that law is fundamen- 
tally altered in a turbulent flow and poses no necessary 
obstacle to fast reconnection. Indeed, there are very close 
relations between the work of Eyink and LV99. Fo r ex- 
ample, it was pointed out in Eyink fc Aluie (2006) that 
the "stochastic wandering" of tield-lines in a rough mag- 
netic field which was invoked in LV99 is exactly analo- 
gous to the "spontaneous stochasticity" phenomenon for 
Lagrangian particle trajectories guided by a rough, tur- 
bulent velocity field. 

In this work we shall establish even deeper connections 
between the work of Eyink and LV99. In particular, we 
show that the detailed, quantitative predictions of the 
LV99 theory can be rederived by applying the stochas- 
tic flux-freezing laws of [Eyink | (|2011| to the problem 
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of magnetic reconnection. In this derivation, the LV99 
theory emerges as the natural turbulent analogue of the 
Sweet-Parker laminar solution, with resistive diffusion of 
field-lines through the reconnection layer replaced by tur- 
bulent Richardson diffusion. To estimate the latter, we 
employ the sa me phenomenological t urbul ence model as 
did LV99, the |Goldreich fc Sridhar"] ( |1995D theory [here- 
after, GS95], although alternative theories would lead 
to similar results. The new derivation demonstrates the 
dynamical consistency of the LV99 theory and also elim- 
inates the need for some difficult arguments developed 
in LV99 concerning the "motion" of already-reconnected 
field-lines. 

A second main purpose of this paper is to explain in 
a clear physical manner the concept of "spontaneous 
stochasticity" , both of Lagrangian trajectories and of 
magnetic field-lines, and the related notion of stochas- 
tic fiux-freezing. "Spontaneous stochasticity" has some 
rather dramatic implications, such as the breakdown of 
Laplacian determinism for classical dynamics. Neverthe- 
less, the phenomenon is not difficult to understand in 
intuitive terms. We shall also try to explain the concept 
of stochastic ffux-freezing in high (kinetic and magnetic) 
Reynolds-number MHD turbulence in the simplest pos- 
sible ma nner, explicatin g some points that were left im- 



plicit in Eyink (2011). An essential point which has 



been emphasized many times be fore ( Newcomb ||1958 
V asyliuna s ||19 7 2[ [A lfven ||1976[) is that the notion o: 
ffux-freezing is largely conventional. While it is a very 
useful intuitive device, the frozenness of ff eld-lines is a 
"meta-physical" principle which cannot be subjected to 
direct experimental test. Any deffnition of ff eld-line mo- 
tion that is consistent with the hydromagnetic equations 
is equally acceptable. We believe that the stochastic ap- 
proach to ffux-freezing will provide a powerful heuristic 
tool in astrophysics. In particular, there should be im- 
portant applications whenever the plasma ffuid velocity 
ffeld has an extended inertial-range and the stochastic 
motion of ffeld-lines is therefore due to turbulent advec- 
tion rather than to resistivity. This paper presents an 
illustrative "case study" of an application to the prob- 
lem of large-scale turbulent reconnection, but m any other 

uses can be envisage d. A s anot her example, |Eyink fc 

Neto I ( |2Q1QD , [Eyink | ( |2Q1Q[ |2QllD discuss the critical role 
ot stochastic tlux-freezmg to the mechanism of the tur- 
bulent kinematic dynamo for magnetic Prandtl numbers 
either small or unity. 

The detailed contents of this paper are as follows: In 
§2 we brieffy recall the just iff cat ion for an MHD descrip- 
tion of astrophysical plasmas. In §3 we discuss the scal- 
ing laws of MHD turbulence and Lagrangian particle dy- 
namics in MHD turbulence, and introduce the concepts 
of stochastic ffux freezing fc spontaneous stochasticity. 
We next discuss the generalized Ohm's law in §4. Us- 
ing the concept of spontaneous stochasticity, we rederive 
the predictions of LV99 theory of turbulent reconnection 
in §5 and consider various limiting cases. We provide a 
discussion of our results in §6 and a ffnal summary in §7. 

2. VALIDITY OF A MAGNETOHYDRODYNAMIC 
DESCRIPTION 

An MHD (or nearly MHD) description of a plasma can 
be justiffed based either on collisionality or on strong 
magnetization. This issue has been well-studied and the 



subject of many discussions, e.g. Kulsrud (1983), so we 
shall here just brieffy review the standard understand- 
ing. There are three characteristic length-scales of im- 
portance: the ion gyroradius pi, the ion mean- free-path 
length irnfp,i arising from Coulomb collisions, and the 
scale L of large-scale variation of magnetic and velocity 
ffelds. Astrophysical plasmas are in many cases "strongly 
collisional" in the sense that imfp,i Pi, e.g. the inte- 
riors of stars and accretion disks. In such cases, a stan- 
dard Chapman-Enskog expansion yields a ffuid descrip- 
tion of the plasma (Braginsky 1965), either a two-ffuid 
model for scales between irnfp,i and the ion skin-depth 
= Pi/ a/A oi" simple MHD at scales much larger than 
Si. However, astrophysical plasmas are often at most 
only "weakly collisional" with irnfp,i ^ Pi, especially 
when they are very hot and diffuse. This can be seen 
from the simple relation 



(X 

Pi 



InA c 



(1) 



which follows from the standard fo rmula for the Coulomb 
collision frequency (e.g. see Fitzpatrick (2011), 
eq.(1.25)). Here A = 47rnA|) is the plasma parame- 
ter, or number of particles within the Debye screening 
sphere. When astrophysical plasmas are very weakly 
coupled (hot and rariffed), then A is extremely large. 
Of course, va < c d^s well, but the ratio va/c is gener- 
ally only moderately small in well-magnetized plasmas. 
Typical values for some weakly coupled cases are shown 
in the following table: 

Table 1 

Representative Parameters for Some Weakly- Coupled Astrophysical Plasmas 



Parameter 


warm ionized 


post-CME 


solar wind at 




ISM" 


current sheets^ 


magnetosphere*^ 


density n, cm~^ 


.5 


7 X 10^ 


10 


temperature T, eV 


.7 


103 


10 


plasma parameter A 


4 X 10^ 


2 X 10^0 


5 X 10^0 


ion thermal velocity Vth,ii cm/ s 


10^ 


3 X 10^ 


5 X 10^ 


ion mean- free-path irnfp,i, cm 


6 X 10^^ 


10^0 


7 X 10^2 


magnetic diffusivity A, cm'^/s 


10^ 


8 X 10^ 


6 X 10^ 


magnetic field B, G 


10-6 


1 


10-4 


plasma beta /5 


14 


3 


1 


Alfven speed va, cm/s 


3 X 10^ 


3 X 10^ 


7 X 10^ 


ion gyroradius pi, cm 


10^ 


3 X 10^ 


6 X 10^ 


large-scale velocity U, cm/s 


10^ 


4 X 10^ 


5 X 10^ 


large length scale L, cm 


1020 


5 X 10^0 


10^ 


Lundquist number Sl = 


3 X lO^s 


2 X 10^^ 


109 


resistive length* i:^, cm 


5 X 10^ 


1 


20 



Norman & ] 



Ferrara |(l996' ; Ferriere )200i; ^ Bemporad et al. |2008; ' Zimbardo et al. j2010 



*This nominal resistive sca le is calculated from ~ L{va/U)S j^'^^ , assuming GS95 turbulence holds 
down to that scale (see eq. |72j ), and should not be taken literally when £^ < Pi- 

The primary expansion to justify a hydrodynamic de- 
scription of these weakly coupled but well-magnetized 
plasmas is an expansion in the small ion gyroradius. This 
yields at all length scales larger than pi a description by 
"kinetic MHD equations", nearly identical to standard 
MHD. The most significant difference is that the pres- 
sure tensor in the momentum equation is anisotropic, 
with the two components p\\ and p± of the pressure par- 
allel and perpendicular to the local magnetic field di- 
recti on deterniined from a 1-dimensional kinetic equa- 
tion (iKulsrud "1983). The magnetic field solves an ideal 
induction equation, if one ignores all collisional effects. 
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although electron-ion collisions produce in reality a non- 
zero resistivity rj. However, in very many cases (cf. the 
ISM and the magnetosphere in Table 1) the resistive 
length-scale is much smaller than pi and also than 

Pe ^ 'is Pi' Magnetic field-lines are, at least formally, 
well "frozen-in" on these scales. Plasmas that are not 
strongly collisional further divide into two cases: "colli- 
sionless" plasmas for which imfp,i ^ the largest scales 
of interest, and "weakly collisional" plasmas for which 
L ^ ^mfp,i' In the latter case the "kinetic MHD" de- 
scription can be further reduced in complexity at scales 
greater than irnfp,i by including the Coulomb collision 
operator and making again a Chapman-Enskog expan- 
sion. This reproduces a fully hydrodynamic MHD de- 
scription at those scales, with anisotropic transport be- 
havior associated to the well- magnetized limit. Among 
our examples in Table 1 above, the warm ionized ISM 
is "weakly collisional", while post-CME current sheets 
and the solar wind impinging on the magnetosphere are 
(nearly) "collisionless." 

Additional important simplifications occur if the fol- 
lowing assumptions are satisfied: turbulent fiuctuations 
are small compared to the mean magnetic field, have 
length-scales parallel to the mean field much larger than 
perpendicular length-scales, and have frequencies low 
compared to the ion cyclotron frequency. These are 
standard assumptions of the Goldreich-Sridhar (GS95) 
theory of MHD turbulence, discussed in the following 
section, but are often valid much more generally in as- 
trophysical plasmas. They are the basis of the "gy- 
rokinetic approximation" which was developed for fu- 
sion plasmas but more recently extensively applied in 
astrophysi cs ( Schekochihin et al. 2007; Schekochihin et 



al. ||2009 ). At length-scales greater than the ion gyro- 
radius Pi, which mainly concern us in this work, a re- 
markable reduction occurs. The incompressible shear- 
Alfven wave modes have an autonomous dynamics un- 
affected by the compressive modes and described by the 
simple "reduced MHD" (RMHD) equations. This fact 
is of fundamental importance for the theory developed 
in the present work, permitting us to base our analysis 
on an incompressible MHD fiuid model. Compressible 
fast magnetosonic waves are subject to strong damping 
both from collisionless effects and in shocks, and, also, 
are not regenerated by the incompressible dynamics (see 
below). The other compressible plasma modes (slow and 
entropy waves) are "passively" transported by the shear- 
Alfven waves. For "collisionless" plasmas the compress- 
ible modes are described by gyrokinetic equations for a 
distribution function g in phase space and dissipated by 
collisionless damping. For "weakly collisional" plasmas 
at scales greater than irnfp,i the slow and entropy waves 
are described by the compressible MHD equations and 
damped primarily by (field-parallel) viscosity and ther- 
mal diffusion. 

Although we are generally interested in this work in re- 
connection phenomena at length scales much larger than 
Pi^ there are some situations where the scales of interest 
are comparable to the ion gyroradius, e.g. in the mag- 
netosphere (see Table 1). At scales smaller than pi but 
larger than p^^ the plasma is described by an ion kinetic 
equation and a system of "electron reduced MHD" (ER- 



et al.||2007| ISchekochihin et al. [|2009D . This system ex- 
hibits the Hall effect", with distmct ion and electron 
mean flow velocities and magnetic fleld-lines frozen-in 
to the electron fluid. The ERMHD equations (or the 
more general "electron MHD" or EMHD equations) pro- 
duce the typical features of "Hall reconnection" such 
as quadrupolar magnetic fle lds in the reconnection zone 
(Uzdensky & Kulsrud 2006[ )p] At length scales smaller 
than kinetic equations are required to describe both 
the ions and the electrons. It is at these scales that the 
magnetic flux flnally "unfreezes" from the electron fluid, 
due to effects such as Ohmic resistivity, electron inertia, 
flnite electron gyroradius, etc. However, as we shall dis- 
cuss at length in this work, these weak effects are vastly 
accelerated by turbulent advection and manifested, in 
surprising ways, at far larger length scales. 

We end this short review with just a few additional ref- 
erences. The effects of compressible modes have been ex- 
tensively studied in n umerical simulations of isothe rmal 
compressible MHD by|Cho fc Lazarianl (|2002[ |2003|) (see 



also Kowal & Lazarian ( |2010|) ). Separatmg MHD fluc- 



t nations into Altven, slow and fast modes, they showed 
that the energy transfer from Alfvenic to compressible 
modes scales oc {Sv)'^ /{v\ + c^), where Sv is the velocity 
perturbation at the given scale and va and Cg are Alfven 
and sound speeds, respectively. As a result, for many 
astrophysically important cases the energy exchange is 
suppressed to a high degree. The Alfven-wave cascade 
then decouples from the compressible degrees of freedom 
and the shear- Alfven modes in the compressible regime 
exhibit scalings and anisotropy very similar to those in 
the incompressible regime. These results give further 
justiflcation in our present study to adopting a simple 
model of an incompressible MHD plasma. As we shall 
see, it is the shear- Alfven modes that are respons ible for 
all the pr i ncipal effects discussed below. Finally, |Cho fc:| 
Lazarian (2004) have carried out a numerical study of 
three-dimensional EMHD turbulence. They observed a 
critically-balanced, anisotropic EMHD cascade quantita- 
tively different, but qualitatively the same, as the GS95 
strong MHD turbulence discussed in the following sec- 
tion. Although we focus in this work on reconnection 
phenomena at scales much larger than p^, similar ideas 
may apply at scales less than pi in the EMHD regime. 

Because the Hall MHD equations have played a prominent 
role in magnetic reconnection research of the past decade (Shayl 



"et al."l998 1999 Wang et al."2000 Birn et al. ,2001, .Draken^UUlf 
Malakit et al. 2009 Cassak et al. 2010), it is wortn remarking 
that those equations are essentially never applicable in astrophys- 
ical environments. A derivation of Hall MHD based on collisional- 
ity requires that the ion skin-depth Si must satisfy the conditions 
6i ^ L ^ irnfp,i- The second inequality is needed so that a two- 
fluid description is valid at the scales L of interest, while the flrst 
inequality is needed so that the Hall term remains sig nificant at 
those scales. However, substituting Si = pi/^/pi into (m) yields the 
result 



A vth,', 
In A c 



MHD) equations for kinetic Alfven waves (| Schekochihin 



The ratio Vth,i/c is generally small in astrophysical plasmas, but 
the plasma parameter A is usually large by even much, much more 
(see Table 1). Thus, it is usually the case that imfp,i ^ ^i, unless 
the ion temperature is extremely low. A collisionless derivation 
of Hall MHD fro m gvrokinetics requires al so a restrictive condi- 
tion of cold ions ( [Schekochihin et al. | ( |200 9), Appendix E). Thus, 
Hall MHD is literally valid only tor cold, dense plasmas like those 
produced in some laboratory experiments, such as the MRX recon- 
nection experiment (Yamada^l999 Yamada et al._2010J . 
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3. MAGNETOHYDRODYNAMIC TURBULENCE 

Magnet ized tu rbulence is a tough and complex prob- 
lem (see Biskamp (2003) and references therein). A 
broad outlook on the astrophysical impli cations of the 
turbu l ence c an be found in a review by [Elmegreen 



Scalol (2004) 
ular clouds 



while the effects of turbulence on m olec- 
and star formation are reviewed in IMc- 



Kee fc O striker (2007) and Ballesteros-Paredes et al. 
(|2007|). However, the issues of turbulence spectrum and 
its anisotropics, we feel, are frequently given less atten- 
tion than they deserve. 

In spite of its complexit y, the turbulent cascade is self- 
similar and universal (see 'Monin & Yaglom (1975)) over 
its inert ial range. The physical variables are proportional 
to simple powers of the eddy sizes over a large range of 
sizes, leading to scaling laws expressing the dependence 
of certain non-dimensional combinations of physical vari- 
ables on the eddy size. Robust scaling relations can pre- 
dict turbulent properties on the whole range of scales, 
including those that no large-scale numerical simulation 
can hope to resolve. These scaling relations are extremely 
important for obtaining an insight into processes on the 
small scales. In the interstellar medium, the inertial 
range of fluctuations is very appreciable, i.e. from hun- 
dreds of k ilometers to dozens of parsecs (see Armstrong 
et^al.|(|1995) for the data at small scale, and Chepurnov & 
Eazarian (2010) for its extension to pc scales) and there^ 
tore the details of the large-scale driving and of the small- 
scale dissipation become unimportant for the properties 
of turbulence at the long range of intermediate scales. 

The presence of a magneti c field makes MHD tur- 
bulence anisotropic (Montgomer y fc Turner| Jl9 81|); 
Matthaeus et al.|(IT983); Shebalin et al.l (I1983D [Hi 



p^eich fc SrldW I ( 11995) , 
(2003) lor a review). 'I'he relative m 



igdc 



see |UughtQn et al 

(2003) for a review). 'I'he relative importance of hydro- 
dynamic and magnetic forces changes with scale, so the 
anisotropy of MHD turbulence does too. If the turbu- 
lence spectrum and its anisotropy are known, then many 
astrophysical results, e.g. the dynamics of dust, scatter- 
ing and acceleration of energetic particles, thermal con- 
duction, can be obtained. 

3.1. The Goldreich-Sridhar picture of strong and weak 
MHD turbulence 

The standard theory for Alfvenic turbulence i s cur- 
rently the one sugg ested by |Goldreich fc Sridhar I (|1995|) 
[henceforth GS95]. The cornerstone ol the theory is the 
so-called critical balance, which is the rough equality be- 
tween the timescale of the eddy-turnovers perpendicular 
to the local magnetic field and the periods of waves as- 
sociated with the eddie^ The predictions of the GS95 
model are in rough agreement with numerical simulations 



(Cho fc Vishniac 2000; Maron fc Goldreich 2001 Cho et 



^ This equality is frequently expressed following GS95 convention 
Tn terms of parallel and perpendicular wave numbers. However, as 
t!ie waven umbers are defined in the global magnetic field frame oT 
referenceTthis may be misleading. Critical balance is a condition 



satisfied only in the frame related to the local magnetic held. The 
original GS95 work also uses closure relations which are true only in 
the global frame of reference. The importance of the local system of 
reference was implicit in Kraichnan ( 1965 ) but was hrst discussed 
withm the GlSyb model m subsequent publications (LVyy, Uho ^ 
Vishniac (2UUU); Mar on «^ Goldreichi(|2UUli);iUho et al.|(|2UU2i);iUKo 
Lazarian (2UU2)j. 



al. 2002 Beresnyak fc Lazarian 2006), although some dis- 
agreement in terms ol the measured spectral slope was 
noted. This disagreement produced a flow of papers with 
suggestions to improve the GS95 model by includ ing ad- 
ditional effects like dynamical alignment (Boldyre v |2005l 



200 6), po larization intermittency (Beresnyak fc Lazar 

|]2006| , non -locality ([Go goberidze 2007). ' " 
cent studies in | Beresnyak fc Lazarian^ ( ,2009( 2010) and 



More 



Beresnyak ( 20 lip indicate that all nunierical simulations 
performed to date may not have sufficiently extended 
inertial range to get the actual spectral slope and there- 
fore worries about the "inconsistency" of the GS95 model 
may be premature. 

We shall add parenthetically that in a number of 
applications the empirical so-called composit e 2D /slab 
mo del of magnetic fluctuations is used (see IBieber et 



al 



(1994)). In the latter model, which is also known 
sis' two- component models it is assumed that interplane- 
tary fluctuations can be described as a superposition of 
fluctuations with wave vectors parallel to the ambient 
large-scale magnetic field (so-called slab modes) and per- 
pendicular to the mean field (so-called two-dimensional 
modes). This model results in a Maltese cross structure 
of magnetic correlations. This model was developed to 
account for the solar wind observations, which it does 
well by adjusti ng the intensity of two components (see, 
e.g., Matthaeus et al. ( 1990| )). From the theoretical point 
of view, as well as from the point of numerical test- 
ing, this standard model is rather vulnerable. Indeed, 
2D fluctuations are consi stent with the theory o f wea k 
Alfvenic turbulence ( see iNg fc Bhattachariee~| (|1996 ); 
Lazarian fc Vishniac| (l999p ; [Galtier et al. [ ( pOOOP j, but 



this analytical theory predicts strengthening ol the cas- 
cade with decrease of the turbulence scale. Therefore 2D 
fluctuations can describe Alfvenic turbulence only over 
a limited range of scales. In addition, slab modes do 
not arise naturally in MHD numerical simulations wit h 
large-scale driving (see Cho fc Lazarian (2002 2003| )). 
We do not discuss this model ol MHD turbulence lur- 
ther, as it does not have physical motivation and is not 
supported by numerical modeling. It may be treated as a 
parameterization of a particular type of magnetic pertur- 
bation dominated by the peculiarities of driving. On the 
contrary, we consider turbulence which has an extended 
universal inertial range over which the effects of driving 
and dissipation are negligible. 

Below we shall adopt the GS95 model of turbulence 
to describe the Alfvenic part of MHD turbulent fluctu- 
ations. Note that the Alfvenic perturbations are most 
important for magnetic field wandering which, accord- 
ing to LV99, is the process that enables fast reconnec- 
tion. The GS95 model can be generalized to compress- 
ible turbulence and even for supersonic motions numeri- 
cal calculations sho w that the Alfvenic perturbations ex- 

Lazarian|20Q2l[2003 



hibit GS95 scal ing (Cho 



Kowal i 



Lazarian||2010| ). In what follows we consider IVlHD tur 
bulence where the flows of energy in the opposite direc- 
tions are balanced. When this is not true, i.e. when the 
turbulence has non-zero cross-helicity, the properties of 
turbulence depart substantially from the GS95 model 

^ Among the existi n g theories of imbalanced turbulence fsee 
[Lithwick et a l.' f2007); Beresny ak fc Laza rian^ ('2008'); 'Chandran 
( |2UU^| ); perez «fe Boldyrev^ ( ^:^UU9| , all, but ^Beresnyak fc Lazarian 
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Similarly, we shall not discuss MHD turbulence at high 
magnetic Prandtl numbers, when the field-parallel vis- 
cosity is much larger than resi stivity (see Cho et al. (2002 
2003| ; Lazarian et al.| (|2004|) . This viscosity-dominatec . 
regime has important consequences for turbulence in par- 
tially ionized gas (Lazarian et al. 2004), while our present 
study deals mostly with conducting fluids and fully ion- 
ized plasmas. Our considerations do apply for large mag- 
netic Prandtl numbers at scales greater than the viscous 
length, if the kinetic Reynolds number is also large. 

The nature of Alfvenic cascade is expressed through 
the critical balance condition in the GS95 model of strong 
turbulence, namely. 



(2) 



where du£ is the eddy velocity, while i\\ and £_\_ are, 
respectively, eddy scales parallel and perpendicular to 
the local direction of magnetic field. The critical bal- 
ance condition states that the parallel size of an eddy 
is determined by the distance Alfvenic perturbation can 
propagate during an eddy turnover. The qualifier "lo- 
cal" is important as no universal relations exist if ed- 
dies are tre ated in respect to the global mean magnetic 
field (LV 99, [Chofc VishniacI (| 2000 ) ; Maron & Goldreich 

0;|Lithwick Goldreich (2001,); |Cho et al.| ( ^2002| )). 
inmg ([2| with the Kolmogorov cascade notion, i.e. 
that the energy transfer rate is Suj/{i±/5u£) = const 



one gets Su£ ~ £]^^ and i\\ ~ i]^^ . If the turbulence 

1/3 2/3 

injection scale is L^, then ex L-^ , which shows 
that the eddies get very anisotropic for small £_\_ . Recent 
measurements of anisotropy in the solar wind are consis- 



tent with th ese predictions ( Podesta ||2010 Wicks et al 
] [2QTq1[2MT] ). 

Gritical balance is a feature of strong turbulence, which 
is the case when the turbulence is injected isotropically 
with velocity amplitude ul = va- If the turbulence is 
injected at velocities ul <^ va (or anisotropically with 
II ^ then the cascade is weak, with i± of the ed- 

dies decreasing but i\\ = Li unchanged by resonant wave 
interactions ([Montgo mery & Matthaeus 



fc Vishniac|l999(|Galtier et al. 



Lazarian 



a result of the weak cascade the eddies get thinner, but 
preserve the same length along the local magnetic field. 
The energ y cascade rate in the weak MHD turbulence 



regime is dKraichnan 1 196 5 Ng & Bhattacharjee 



Lazarian Vishniacll 



1996 



r5uj/il^SujLi/vAil, 



(3) 



where Tcorr = Li/vA is the decorrelation time due to 
large-scale Alfven waves. This implies that the scaling of 
weak turbulence is 

(4) 



([2 008 ) seem to contradict to numerical testing in 'Beresnyak &' 
^izarian (2009 2010). Solar wind presents a system with imbal- 
anced turbulence. In compressible media the imbalance decreases 
due to reflecting of waves from pre-existing density fluctuations 
and due to the development of parametric instabililites. 

^ To stress the difference between local and global systems here 
we do not use the language of fc- vectors. Wavevectors parallel and 
perpendicular to magnetic flelds can be used, if only the wavevec- 
tors are understood in terms of a wavelet transform deflned with 
the local reference system rather than ordinary Fourier transform 
deflned with the mean fleld system. 



For weak turbulence initially the LHS of Eq. ([2| is larger 
than the RHS, but as decreases this eventually makes 
Eq. ([2| satisfied. 

Comparing Eqs. (|4| and ([2| for = L^, one can see 



that the transition to the strong MHD turbulence hap 
pens at the scale L trans = LiM\ and the velocity at 
this scale is v trans = ulMa^ with Ma = ut Jva ^ 1 
the Alfvenic Mach number of the turbulence ([Lazarian 



& Vishniac 1999 Lazarian 2006). Thus, weak turbu- 
lence has a limited, i.e. [Li^LiMj^] inertial interval be- 
fore it gets into the regime of strong turbulence at smaller 
scales. Note that weak and strong are not the character- 
istics of the amplitude of turbulent perturbations, but 
the streng th of non- line ar interactions (see more discus- 
sion m ^Cho et al.|(2003| )) and very weak Alfvenic pertur- 
bations can correspond to a strong Alfvenic cascade. 

While GS95 assumed that the turbulent energy is 
isotropically injected with amplitude ul = va the 
scale L^, LV99 provided general relations for the turbu- 
lent scaling at small scales for the case that the injection 
velocity ul'i^ less or equal to va- Combining Eqs. p|, ([3| 
at = and the constant fiux condition, we get the 
relations between the parallel and perpendicular scales 
of eddies in the strong GS95 cascade range, which can 
be written in terms of £\\ and £± (LV99): 



(5) 
(6) 



As we discuss later, the present day debates of whether 
GS95 approach should be augmented by additional con- 
cepts like "dynamical alig nment" , "polarization" , "non- 



locality" ([B oldyrev 



2009 Gogoberidze 2007r 



Beresnyak & Lazarian 2006 
o not change the nature of the 
reconnection of the weakly turbulent magnetic field. 

The above relations were exploited in LV99 to esti- 
mate magnetic field lateral diffusion. Moving a distance 
s along field-lines, one finds that a pair of lines initially 

a distance apart at 5 = separate at the rate 

d 5be Sue 
ds Bq Va 

Substituting the scaling given by Eq. (|6|, one can solve 
to obtain 

.2 .W.3/r.AA^4 (g) 



when Li > s ^ . This result applies only until s = Li 
when i±{s) Ltrans-, the transitional scale between weak 
and strong turbulence. Since l\\ = Li is constant in weak 
turbulence, there are no eddies with parallel length-scale 
greater than Li. Thus, as s continues to increase, the lines 
encounter at every increment of 5 by a new turbulent 
eddy and undergo a further separation by Ltrans- The 
result is a diffusive random-walk with 

4 ^ Llansis/Li) ^ sLiM\ (9) 

for s > Li. The results (|8|,(|9| are exact analogues 
for magnetic field-lines of the 2-particle diffusion of La- 
grangian particle trajectories discussed in the next sec- 
tion and can be obtained more rigorously by the methods 
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employed there. Since this discussion is somewhat out of 
logical order, we defer it to Appendix A. 
A remarkable feature of ([8| , Q is their complete lack of 

dependence on the initial separation of the lines. As 
we shall see below, it is this property which leads in the 
LV99 model to fast magnetic reconnection in MHD tur- 
bulence independent of the microscopic resistivity. The 
shear- Alfvenic component of turbulence i s the most im- 
por t ant for mag netic field-line wandering. |Cho fc: Lazar- 



lan 



( 2002 , 2003 ) demonstrated the correspondence of the 
Alfvenic component of compressible MHD turbulence, 
obtained by decomposition into fundamental modes, to 
the Alfvenic turbulence in incompressible media. Thus 
we shall use the relations above to treat magnetic fields 
in realistically compressible astrophysical fluids. 

3.2. Lagrangian particle dynamics in MHD turbulence 

It will be important for what follows to develop a theo- 
retical understanding of the dynamics of Lagrangian par- 
ticles in MHD turbulence. Here one means by "particles" 
not the microscopic charged constituents of the plasma 
(electrons and ions), but instead macroscopically small 
and microscopically large parcels of plasma fluid. The 
Lagrangian perspective is particularly important to un- 
derstand turbulent processes and there has been an ex- 
plosion of resear ch in ^t his are a over the last deca de (see 
Falkovich et al. (|200 1 1 ; [Salazar fc Collins | ( |2009| ) ; [Toschi 

|). 'inhere has been relatively little 

work done on the Lagrangian properties of MHD turbu- 
lence, however, where the effects of Alfven waves on fluid 
particle motion can be signiflcant. We shall here develop 
a theory of Lagrangian MHD turbulence at the level of 
GS95 phenomenology by a semi-quantitative use of the 
equations of motion. 

We begin with the simplest case of 1-particle or Taylor 
diffusion, which concerns the displacement Jx(t) = x(t) — 
xo of a Lagrangian fluid particle from its initial position 
as it moves according to the advection equation 



dt 



x(t) = U(t) 



(10) 



with U(t) = u(x(t),t) the Lagrangian particle veloc- 
ity. At early times t <C or short compared with the 
eddy-turnover time tl = ul/Li^ m otion is ballistic with 
{x'^{t)) - u\t^. It was shown by [Taylor | (|l92l|) that 
the motion is diffusive {5xi{t)5xj{z)Y 
t^TL^ with the eddy-diffusivity 



2Vijt at times 



+ 00 



dt{Ui{t)Uj{0)). 



(11) 



These results carry over unchanged to MHD turbulence. 
For the weak magnetic fleld case with Ma > 1 (including 
hydrodynamic turbulence with Ma = cxo), {U{t)U{0)) 



-t/TL 



Thus, 



( [Mordant et al. ||200l[ [Busse fc Muller [|2008[ ). 



D 



(12) 



exactly as for hydrodynamics of neutral fluids. When 
Ma — 1, one expects slight differences from the hydro- 
dynamic case due to anisotropy of MHD turbulence. For 
example, the velocity components U\\ and ought to 
have times scales tl\\ > tl^^ since the parallel compo- 
nents are pseudo- Alfvenic modes with essentially passive 



scalar dynamics, due to reduced nonlinear interactions. 
Thus, one exp ects that D\\ > D± : as indeed observed in 
simulations of "Busse fc Miiller I (|2008|). For <C 1, 



there are more profound efl:'ects due to the strong fleld. 
Weak MHD turbulence becomes relevant, for which the 
nonlinear time-scale is tl = M^^uj^ with uja = vaI^i 



the forcing-scale Alfven wave frequency. (See K raichnan 
II (119651) ; iSridhar fc Goldreichj ([1994|); [Galtier et al. 
(|2(]0nin ct. also eq.(|25| below foFTT = Li), l^urther- 
more, fluid particles m the Lagrangian frame experience 
fast shear- Alfvenic oscillations in velocity, so that 

(f/±(t)f/±(0)> ~ (e'--*-*/--) . (13) 

Thus, ([n]) gives 

(14) 



{uatlY^I 



Since ujatl = M^^ > 1, 



D 



1 



UJaTL 



ulLM%. 



(15) 



There is a large reduction in the turbulent 1-particle 
diffusivity due to inefficient advection by Alfven wave 
modes. Physically, the velocity of a particle transported 
by Alfven wave turbulence experiences rapid oscillations 
in sign that lead to large cancellations in the net dis- 
placement. 

We note as an aside that the Taylor 1-particle diffusiv- 
ity given by (11) coincides with the effective turbulent 
diffusivity of an advected scalar fleld (e.g. temperature), 
whenever the "eddy-diffusivity" concept is valid. In or- 
der for such a description of the turbulence effects to 
be accurate, the scalar flelds must vary slowly on spa- 
tial scales of order Li and, in that case, they experience 
an augmentation of their molecular diffusivity by a "tur- 
bulent diffusivity" precisely give n by ([ll|) - Q^r results 
(p!2|),(15) thus reproduce those of Lazarian (2006) for the 
e3!ay tliermal diffusivity in MHD turbulence. Of course, 
in general, when there is no separation between the scale 
of variation of the scalar fleld and the scale Li of the tur- 
bulence, th en the gradient- transp ort assumption breaks 



down (e.g. Tennekes fc Lu mley [ ([T972), §2.3). In such 



cases, eddy- diffusivity models can produce erroneous re- 
sults and may suflice only for order-of-magnitude esti- 
mates of scalar transport. 

Two-particle turbulent diffusion or Richardson diffu- 
sion concerns instead the separation i(t) = x(t) — x'(t) 
between a pair of L agrangian fluid particles. It was pro- 
posed by Richardson | (1926) that this separation grows 
in turbulent tlow according to the formula 



dt 



(16) 



with a scale-dependent eddy-diffusivity D{i). In hydro- 
dynamic turbule nce Richardson deduced that D{£) ~ 
^1/3^4/3 j^ggg also[Obukhov]([l94l[)) and thus i^{t) - et^. 



An analytical formula for the 2-particle ed dy-diffusivity 
was later derived by ,Batchelor , (1950] ); [Kraichnan 



(1966): 



J — cx 



dt{6Ui{e,o)sUj{e,t)) 



(17) 
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with SUi{i^t) = Ui{'K-\-i^t) — Ui{:si^t) the relative velocity 
at time t of a pair of fluid particles which were at posi- 
tions X and x + £ at time 0. This formula is the analogue 
for 2-particle eddy-diffusivity of Taylor's formula ( 11 ) for 
1-particle diffusivity and it is valid in MHD turmllence 
just as in hydrodynamic turbulence. We can assume sim- 
ilarly as before that 

where ua/ is the Alfven wave- frequency and is the 
nonlinear interaction time both at length-scale i. The 
result for the 2-particle diffusivity analogous to ( 14 ) is 



(19) 



We consider here only the case Ma < 1, for which weak 
turbulence holds when i± > Ltrans and strong turbu- 
lence when i± < Ltrans ' We shall also treat particle dis- 
persion only in the direction perpendicular to the back- 
ground magnetic field, since this is what i s required later 
f or the a pplicatio n to reconnection. See [Busse fc Miiller 
[] (2008); Eyink |2011, ) for some discussion of particle 
dispersion in the field-parallel direction. 
We first consider the strong turbulence regime when 

i±{t) < Ltrans- this casc, T£ = e~^l'^^l^ and a;^^ = 
vaI^W^ where is the parallel scale of the turbulent eddy 
with perpendicular scale ti_(t) (and not the particle- 
separation (t) in the field-parallel direction). Because 
of the critical balance condition in strong MHD turbu- 
lence, 

^Ai^i = {const.) (20) 



In that case 



{ujAinY 



fn 



(21) 



with / = [{(jOaitiY + < 1 a constant factor which 
represents the reduced efficiency of particle transport due 
to wave oscillations. Eq.(19) gives 



(22) 



with the same scaling as the original Richardson | ( |1926[ ) 
2-particle diffusivity. Thus, 



(23) 



showing the same t^-law as in the hydrodynamic case, 
with only a reduced coefficient. 

We now consider the weak turbulence regime with 
^^(t) > Ltrans- The pair of particles are advected apart 
by a turbulent eddy with corresponding perpendicular 
length-scale but with £\\ = Li for all eddies, so that 
00 A £ = = VA/Li^ independent of £±. The nonlinear 
interaction time is estimated from 



n 

together with Q to give 



vaLi 



Since 



n = {i±/vA)Mf. 

^An = Ltrans > 1, 



(24) 

(25) 
(26) 



it follows from ((19]) that 



ULLiM\. 



(27) 



Note that D±{i) is identical with the 1-particle diffusiv- 
ity D± given by ([l5| and, in particular, is independent 
of £± . It thus follows that in the weak turbulence regime 



(28) 



and 2-particle diffusion is identical to 1-particle diffu- 
sion. It should be possible to verify this interesting pre- 
diction by analytical methods of weak turb ulen ce theory 
for Lagrangian particle motion (see B alkl ( |2Q02|). The 
crossover between the dispersion laws and (|28[) oc- 
curs at the time tc ^ Li/vA = MAth- 

4. FLUX-FREEZING AND SPONTANEOUS 
STOCHASTICITY 

There is an intimate connection between magnetic 
reconnection and flux- freezing. Indeed, some authors 
consider that reconnection is precisely any topology- 
changing evolution of the magnetic field due to nonex- 
istence o f a smooth, fiu x-preserving velocity field. For 
example, I Greene (1993|) proposed to define reconnection 
"as evolution m which it is not possible to prese rve the 
global identifi cation of some field lines." See also Hornig 
fc Schindler | (p^96). For this reason, a necessary prelude 
to our treatment of turbulent magnetic reconnection is a 
careful discussion of how fiux-freezing operates in MHD 
turbulence. 

4.1. Stochastic flux-freezing 

Magnetic field lines do not move. Magnetic fields 
evolve under the Maxwell equations, of course, but sin- 
gle field lines do not possess an identity over time. As 
has long been understood, magnetic field-line motion is a 
"metaphysical" concept, defi ned by convention and not 
testab l e by e xperimen t 



Unas 



(1972); Alfven 



leads to agreement wit! 



(cf. Newcomb (1958); Vasyli- 



(1976 )). Any line-velocity which 
1 the Maxwell equations for the 
magnetic field evolution is equally acceptable, e.g. the 
bulk plasma velocity u in ideal MHD, the ExB drift ve- 
locity in low energy plasmas, etc. Flux-freezing is such 
a powerful heuristic tool, however, that it is often for- 
gotten that the concept of magnetic line-motion rests on 
certain conventions and is not a direct, physical real- 
ity. In particular, real physical plasmas always possess 
a Spitzer plasma resistivity (along with other possible 
forms of non-ideality, such as electron inertia, electron 
pressure anisotropy, etc.) The validity of fiux-freezing in 
the standard sense depends upon careful consideration 
of the limit when such non-ideal effects are taken to be 
small. 

Since field-lines do not move in reality, one is free to 
adopt any convention for their motion which is consistent 
with the dynamical equations. Let us for the moment 
assume the validity of the resistive induction equation 



dt'B = Vx(uxB) + AAB 



(29) 



with A = rjc^ /At: the magnetic diffusivity for a simple 
scalar resistivity 77. A smooth, deterministic velocity 

^ It may be that the flux-freezing constraint is broken in actual 
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field is flux-preserving if the magnetic field obeys 

StB = Vx(u*xB). (30) 

Unfortunately, for the resistive MHD model in three 
space dimensions, there is general ly no fiux-preserving 
v elocity whatsoever. For example, |Wilmot- Smith et al. 
U (2005) provide a counterexample with a closed mag- 
netic tield-line through which the flux is changing in time, 
so that no deterministic flux-preserving velocity can ex- 
ist. This seems to leave one without any consistent and 
meaningful approach to discuss field-line "motion" in the 
presence of Ohmic resistivity. 

On the other hand, magnetic line-motion in a resistive 
MHD plasma m ay be very naturally regarded as stochas- 
tic (Eyin k [2009 ). Let the plasma fluid velocity be u(x, t), 
assumed here to be divergence-free for an incompressible 
plasmaj^ Then we may define stochastic Lagrangian tra- 
jectories x(a, t) by the initial- value problem 



f(a,i) 
x(a,io) ■ 



■- u{x{a,t),t) 
a, 



2Xr]{t), 



(31) 



where rj{t) is a three-dimensional vector white-noise pro- 
cess, delta-correlated in time. For any initial magnetic 
field Bo(x) one can form the random field 



B(x,t) = Bo(a).V,X(a,t)|5 



a(x,t) 



(32) 



where a(x, t) is the spatial inverse map to x(a, t). If the 
flow map were deterministic, then (32 ) would be the stan- 
dard Lundquist formula for the magnetic field. What is 
true in the present case with A > is that the ensemble- 
average field 

B(x,t) = B(x,t) (33) 

over all independent realizations of the white-noise is the 
exact solution of the resistive induction equation (29) 
with initial condition 

B(x,to) = Bo(x). (34) 

In fact, (mi-rt33f are mathematically equivalent to ([29| 
and (34)~lt sliould not be too surprising that averag- 



ing over Brownian motions will reproduce the Laplacian 
diflFusion term AAB. 



The eqs.(31 )-(33) are equ ivalent to the following path- 
integral formulaTaerived in Eyink ( |2011 ): 



B(x,t) 



Pa Bo(a(to))-J(a,t) 



a(t)=x 



exp^-^^ |a(r) -u^(a(r),r) 



(35) 



where B is interpreted as a 3-dimensional row vector and 
J7(a, r) is a 3 X 3 matrix satisfying the following ODE 

astrophysical plasmas by some other microscopic mechanism (e.g. 
electron inertia, finite electron gyroradius, etc.). We shall consider 
such alternatives below (section 4.3) and show that they do not 
alter our final conclusions. In particular, we argue that the large- 
scale global reconnection rate in a turbulent plasma is independent 
of whatever microscopic plasma process produces field-line break- 
ing. It is thus pedagogically useful to begin with the simplest case. 

^ The assumption of incompressibility is convenient to simplify 
the arguments and, also, realistic when considering the RMHD 
dvn a mic s of shear- Alfven modes. As a matter of fact, the formulas 
(31|)-(|37| below all generalize to the compressible case. See [ Eyink | 

(wni'r^ 



start 





average 




Fig. 1. — Illustration of the stochastic Lundquist formula. Three 
stochastic Lagrangian trajectories running backward in time from a 
common point are shown in red, blue, green. Starting field vectors, 
represented by correspondingly colored arrows, are transported 
along the trajectories, stretched and rotated, to the common fi- 
nal point. These "virtual field vectors" are then averaged to give 
the resultant magnetic field at that point, indicated by the black 



along the trajectory a(r): 
\ J{a,to)=I. 



(36) 



By applying Va to (|31| it is easy to see that Jia, r) 
Vax(a(io),T) and thus we may identify 



B(a(r),r) = Bo(a(to))-J(a,T). 



(37) 



As often with such path-integration techniques, only the 
initial field Bo(x) and the final field B(x,t) have real 
physical meaning and the intermediate random fields 
B(x, t) are merely calculational devices. We shall re- 
fer to the latter as the "virtual" magnetic fields, since 
they enter as non-real, intermediate states between two 
physical fields at different times. Equation (33) [or (35)] 
shows that summing over all of the "virtuar fields at 
any instant reproduces the true, physical magnetic field. 
It follows from (|32| [or ^] that each "virtual" field 

B is topologically equivalent to the initial field Bq and 
is simply stretched, deformed and advected by its own 
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stochastic flow x. Field-lines for two different "virtual" 
fields that arise from statistically independent flows may 
pass directly through each other as they move, since they 
correspond to completely different random realizations. 
The topology of the field-lines is only changed by the fi- 
nal averaging step, which corresponds to resistive gluing 
of all the "virtual" magnetic field vectors which arrive 
to the same point x at time t. It is this resistive aver- 
aging which reconnects the lines of the initial magnetic 
field Bq. See Fig. 1 for a pictoral representation of this 
process. 



We shall refer to either the relatio ns (|31| )-(33) or the 
equivalent path-integral formula (35)-(|36p as tne stochas- 
tic flux-freezing (SFF) relations. TTiey provide a consis- 
tent approach to describe the "motion" of magnetic field- 
lines in the presence of resistivity. Setting A = one for- 
mally recovers the standard deterministic flux-freezing 
relations of ideal MHD, but this conclusion is not rig- 
orous. A careful discussion is required of the limiting 
process involved. 

4.2. Spontaneous stochasticity 

Standard flux-freezing arguments are invoked in the 
limit of high conductivity or vanishingly small magnetic 
diffusivity. It is usually assumed that flux-freezing in the 
conventional sense will hold better and better as A ^ 0. 
This is not true, however, when the plasma flow is tur- 
bulent. Instead, flux- freezing in this circumstance breaks 
down in a compl etely novel and une xpected way: it re- 
mains stochastic ( [Eyink ||2007| |201lD . 

Consider first the simpler situation of a smooth, lami- 
nar solution of the resistive MHD equations in the limit 
of very small A. Define the standard (deterministic) La- 
grangian flow map by 



f (a,t) = u(x(a,t),t), 
x(a, 0) = a. 



(38) 



It is not hard to see under the stated assumptions thalj^ 
|X(a, t) - x(a, t) I = 0{VXt) (39) 

for typical realizations of the "virtual" flows x. The real 
magnetic field is also given to very good accuracy by the 
standard Lundquist formula, with errors that vanish as 
A ^ 0. Thus, each "virtual" field line is wiggling a small 
distance ~ a/A^ around a particular physical field-line 
for very small A (and not too large times t). For this 
reason, one does not need to make use of the concept of 
"virtual" field lines for smooth, laminar solutions of the 
MHD equations. The lines of force of the real, physi- 
cal magnetic field are themselves "nearly" frozen-in for 
very small resistivities and one does not need to make an 
essential distinction between real and "virtual" fields. 

This is not the case for rough, turbulent solutions of the 
resistive MHD equations, even for vanishingly small A. 
By "rough" fields we mean here more precisely that u, B 

^ A more precise statement is that {|x(a, — x(a, < 
^(e^* — 1) where K is the "Lipschitz constant" of the veloc- 
ity field u. See [Freidlin Went zeh (1984), Ch.2. In the presence 
of Lagrangian cJiaos this upper bound may be attained, with K 
the leading Lyapunov exponent. Note that for t <^ 1/K the upper 
bound reduces to 6At, the usual diffusive estimate. The important 
fact is that the bound vanishes as A — )► 0. 



which solve the MHD equations have power-law energy 
spectra ex with 1 < n < 3, for a long range of 

wavenumbers k with kinematic viscosity u and magnetic 
diffusivity A both small. For example, in the GS95 theory 
of MHD turbulence the three-dimensional (anisotropic) 
energy spectra for both u, B can be written in the form 



i?(fc±,fc||)~t;lM7%: 



4/3^-1/3,-10/3 



/ 



k\\L, 



1/3 



■ (40) 



where we added, compared to the original expressions, 
the dependence on Ma following LV99. The original 
GS95 theory was formulated for Ma = 1, as we men- 
tioned earlier. Note, that the above expression is ex- 
pected to be valid for large enough k\\ and k±. The func- 
tion / is not specified in GS95, but Cho, Lazarian & 
Vishniac (2002) showed that it can be fitted by a Cas- 
taing function (see Castaing, Gagne & Hopfinger 1990), 
which, in its turn, outside the vicinity of /cy = 0, can be 
approximated by an exponent. As we mentioned earlier, 
A: 1 1 and k± in the expression are not the wavevectors in 
the conventional sense as the corresponding scales should 
be measured in the local system of referenc^ 

One thus finds that E{k_\_) ex kj^^^ and E{k\\) cx /c^^, 

so that the fields predicted by GS95 are "rough" in our 
sense. In this case, (39) is invalid except for very short 



times, due to the properties of 2-particle Richardson dif- 
fusion. For example, in the transverse direction perpen- 
dicular to the local magnetic field, (39) holds only for 
times less than tx = 0{^\/£). If one considers a pair of 
independent "virtual" flows x, x', then for times t ^ tx 
in the root-mean-square sense 



|5i(a,i)-51(a,i)|~(£t3)i/2. 



(41) 



See eq.([23|). This result is completely independent of 
magnet icdiffusivity, with very remarkable consequences. 
As X^u 0, and thus also ^ 0, the two realizations 
stay far apart for all times t and do not collapse to a sin- 
gle real particle trajectory. This phenomenon has been 
called spontaneous stochasticity because the distance be- 
tween "virtual" particle trajectories stays finite (and ran- 
dom) in the limit, similar to the way that a spontaneous 
magnetization develops in a ferromagnet below the Curie 
tempe rature when an external magnetic field is taken to 
zero ([B ernard et al. 1998; Gawe dzki & Vergassola 2 0001 
IE vanden Eijnden 2000, 200 1; . Chaves et al. 2003 ; ,gu-l 
piainen 2003). In the limit A, all of the "virtual" 
flows x(a, t) solve the deterministic initial- value problem 
(38), which then possesses infinitely many solutions! For 



more discussion of spontaneous stochasticity and for a 
review of experimental and numerical evidence of this 
remarkable phenomenon, see [Eyink | ( |2011| ). 

It is an immediate consequence that m a "rough" tur- 
bulent velocity field, the conventional notion of flux- 
freezing must break down. The violations are not small. 
Taking the square of (39) as the conventional estimate 
for line-slippage, (^^) frozen ^ At, and comparing with 

^ The use of k vectors in GS95 stems from the fact that the 
concept the local system of reference was not introduced till the 
later works. 
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the predicted amount (41), ~ et^ ^ one finds that 

2 

5l, (42) 



(a) 



frozen 



Ml 



t 



tL 



where II = Li/u^ is the eddy-turnover time and Sl = 
vaLi/X is the Lundquist number based on the turbu- 
lence injection scale. For Sl ^ ^ typical of astrophysi- 
cal systems, the violations of standard flux- freezing are 
enormous and increase with time. It is also not true how- 
ever that flux-freezing is completely violated. Stochastic 
flux- freezing in the sense of (31)- (33) remains valid for 



the "virtual" particle trajectories in a plasma flow which 
corresponds to a high (kinetic and magnetic) Reynolds- 
number, turbulent MHD solution. 

There is another way to understand spontaneous 
stochasticity by considering the real trajectories of ac- 
tual plasma fluid elements, which obey (381). Suppose 
that one considers two fluid-elements located initially at 
locations a, a' which are displaced from each other by a 
small amount p = a^ — a. At early times th e separation of 
the particles is ballistic (Batchelor "1950'). For example, 
GS95 theory implies that in the transverse direction 

|x^(a^t) - xx(a,t)| - |ux(aO - ux(a)|t - {ep^y^H 

(43) 

in the rms sense, for times t < to = O [{p\/e)^^^) and 

separations p± > p^ = 0{{i'^ /e)^^"^). However, for 
the Richardson-type result 



|xx(a^^)-x^(a,^)|-(£t^) 



3x1/2 



(44) 



again holds in the rms sense. Note that the initial parti- 
cle separation p is completely "forgotten" at sufficiently 
long times! See Fig. 2. In the limit taking first v ^ {) 
and then p ^ 0, one obtains infinitely many solutions of 
the deterministic initial- value problem (38). In fact, one 
can see by taking a' = a+p for a randornaisplacement p 
that the real fluid particle trajectories also become ran- 
dom in this limit. More precisely, consider an ensemble 
of random displacements inside the ball [p\ < po- Then 
one obtains a random ensemble of real particle trajecto- 
ries by solving (38) for x(a + p, t) and taking the limits 
first u ^ ana then po 0. In fact, it can be shown 
using Richardson's diffusion approximation for 2-particle 
dispersion that this is the same as the random ensemble 
of "virtual" particle trajectories obtained by solving JsTt 
and taking the limits u^ X ^ together (E & vanden Eijn- 
den |2QQQ| . Put another way, for every "virtual" particle 
trajectory x(a, t) there is a real fluid particle trajectory 
x(a-hp, t) such that |p| ^ and |x(a, t) — x(a+p, t)\ ^ 
in the limit X 0. "Virtual" particle trajectories coin- 
cide with real fluid particle trajectories in this statistical 
sense. 

The breakdown of conventional flux-freezing due to 
"roughness" of the fields and spontaneous stochasticity of 
the Lagrangian flows means that magnetic line-topology 
is no longer preserved in time for the limit A ^ 0. The 
infinite ensemble of "virtual magnetic fields" frozen-in 
to the stochastic flows all have exactly the same line- 
topology as the initial magnetic field. However, the av- 
erage over the ensemble of virtual field-lines that arrive 
to the same final point which gives the resultant physi- 
cal magnetic field will, in general, change the topology. 




(b) 




(c) 




Fig. 2. — Illustration of spontaneous stochasticity. Three pairs 
of Lagrangian trajectories are shown in panels (a),(b),(c), with ini- 
tial separations at time between the two members of the pair 
indicated by arrows and final separations at time t indicated 

by red arrows. The initial separations are progressively decreased 
going from (a) to (b) to (c), but the final separations remain the 
same! Continuing in this manner, one obtains as a limit the situa- 
tion with more than one Lagrangian trajectory for the same initial 
point. 



In particular, this permits fast topology change (mag- 
netic reconnection) independent of the precise value of 
the resistivity. It is worth noting, however, that not all 
topological conservation laws are vitiated. Magnetic he- 
licity conservation, for example, is particularly robust, 
even when velocity and magnetic fields ar e ext r emely 
"rough". It has been proved by Caflisch et al. (1997) 
[Theorem 4.2] that magnetic helicity is conserved m the 
limit as A ^ even for such "rough" fields, as long as the 
scaling exponents of velocity and magnetic field 3rd-order 
structure functions remain positive. This result is con- 
sistent with very extreme singularities such as sparsely 
distributed tangential discontinuities (current and vortex 
sheets). Thus, fast reconnection due to the mechanism 
of spontaneous stochasticity will generally conserve mag- 
netic helicity. 

One last important remark: everything that we have 
stated about Lagrangian particle trajectories applies just 
as well to the magnetic field-lines ^(x, cr), defined by solv- 
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ing (at each fixed time t) 



B(^(x,a),t), 



C(x,0)=x. 



(45) 



Here ^(x, a) is tiie field-line which passes through point 
X (at time t). The parameter a is related to arc- length 
s along the field-line by ds = |B(^(x, a), t)|(icr. In GS95 
theory, field lines which correspond to two nearby points 
X, x' with r = x' — X will separate in the transverse di- 
rection as 



IC±(x',s)-^^(x,s)|~(sVL, 



Nl/2 



Mi 



(46) 



in the rms sense, for s sq = 0{vA{r'j_/e)^^^) and for 
> rx = O {{X^/ey^^) . See eqM. This is exactly 
the "stochastic line- wandering" whicn was invoked in the 
LV99 theory of fast magnetic reconnection and, later, in 
the theory of thermal conduction in a tur bulent MHD 
plasma (see Narayan & Medvedev"! (|2001') for Ma = 1 
and Lazarian (2006) for an arbitrary Ma)- Just as for 
Lagrangian trajectories, the initial separation r between 
field lines is "forgotten" after proceeding a large enough 
distance along them. This is the reason that the afore- 
mentioned results for the thermal conductivity do not 
depend upon the electron Larmor radius pe- In fact, mag- 
netic field-lines become stochastic in the limit A — ?► in 
precisely the same sense as do Lagrangian particle trajec- 
tories in the limit u oF^ Consider a random bundle of 
field lines ^(x + r, s) for a stochastic displacement vector 
r distributed over the ball |r| < tq of radius tq. By taking 
the limits first A ^ and then ro ^ 0, one obtains an in- 
finite ensemble of magnetic field-lines all passing through 
the same point x! 

4.3. The generalized Ohm^s law 

Our arguments for the essential stochasticity of tur- 
bulent fiux- freezing in sections 4.1-2 may have given the 
impression that the origin of the randomness is necessar- 
ily connected with Ohmic resistivity. If our conclusions 
really did depend upon the resistive MHD model (29) as- 



sumed there, then they could be criticized as physically 
unrealistic. Indeed, there are in actuality many non-ideal 
terms that appear in the magnetic equations of motion 
for an ionized plasma, which may be summarized in the 
generalized Ohm's law 



E = -uxB + T^^J^ 
nec ne^ V dt 



JxB 



nec 



V-(uJ + Ju 



JJ) (47) 



when quasineutr ality and <C m ?; are as s umed. See | Va- 
syhunas ( 1975]); [Priest fc Forbes | ( |2000[ ); [ Bhattacharjee 



et al. I (|1999| ) .T'he electric tields appearmg on the right- 



hand side are, respectively, the motional field, Ohmic 
fields associated to perpendicular and parallel resistivi- 
ties, the Hall field, a contribution from the electron pres- 
sure tensor, and electron inertial contributions. All of 
these terms have been invoked in various theories of fast 
magnetic reconnection. However, we shall argue that 

In the limit of large viscosities that we do not consider here 
the physics of field wandering is somewhat different and requires a 
separate discussion (see Lazarian, Vishniac & Cho 2004). 



in the turbulent environments that we consider none of 
these terms alters our fundamental conclusions. In the 
first place, the stochasticity of flux- freezing in turbulent 
plasmas is due not to resistivity (or to other terms in the 
generalized Ohm's law) but instead to the "roughness" 
of the MHD solutions in a long inertial range. In the sec- 
ond place, the precise microscopic plasma mechanism of 
"line-breaking" that acts at small scales (below the ion 
or electron gyroradii) is irrelevant to the inertial-range 
turbulence dynamics, which will be fundamentally the 
same for any such mechanism. 

To demonstrate the first point, we show that flux- 
freezing is effectively stochastic in the turbulent inertial- 
range even for the ideal induction equation 



StB = Vx(uxB) 



(48) 



with all non-ideal terms in ( |47[ ) set to zero! Here we as- 
sume that u, B are smooth at length-scales below some 
cut-off < £d but rough (turbulent) at larger scales. Notice 
that the above equation (|48| is realistic at leading-order 
for collisionless, magn etizea plasmas at scales larger than 
the ion-gyroradius pi ( Kulsrud |1983| , while the smooth- 
ness scale id set by field-perpendicular viscosity and re- 
sistivity is often of the same order or smaller than pi. 
Thus, our assumptions are quite realistic. Because of 
the cut-off id, flux-freezing in the standard sense must, 
in fact, be valid. How then can we claim that it becomes 
stochastic? To see this, consider the magnetic field ob- 
served at some finite space-resolution p : 



Bp(x,t) = / d^r Gp(r)B(x + r,t), 



(49) 



where we have introduced a coarse-graining kernel Gp 
to represent the smearing effect of the observation over 
a ball of radius p around the space point x. We shall 
assume below that id <^ p <^ Li, the scale of the largest 
turbulent eddies. Thus, p^ pi satisfies these conditions. 
Applying the standard Lundquist formula for the frozen- 
in magnetic field, one obtains 



Bp(x, 



B(a,to)-VaX^,^o(a) 



det(VaXt,to(a)) 



Xt,tQ (a)=x+r 



Gp{v).d'r 
(50) 

The Lagrangian particle trajectories that appear in this 
formula start in the ball of radius p around x at time t 
and then follow the flow velocity u backward in time to 
to, as illustrated in Fig. 3. When p lies in a GS95 turbu- 
lent inertial range, then the trajectories explosively sepa- 
rate to a perpendicular distance Ax^ ~ (£|t— top) in- 
dependent of p at times |t — to| ^ (p^/e)^/^. The result is 
indistinguishable from the stochastic Lundquist formula 
(35) which was derived in section 4.1 using the stochastic 



representation of Laplacian resistivity. In fact, in a for- 
mal mathematical limit takin g fir st id 0, then p ^ 0, 
the Lagrangian trajectories in (50) remain stochastic and 



the two formulas coincide. This is a rigorous theorem for 
the Kazantsev-Kraichnan dynamo model, where it has 
been proved that the ensemble of stochastic Lagrangian 
trajectories as constructed above is precisely the same as 
that obtained for the A ^ limit (E & vanden Eijnden 
|| |2QQQ| ). Stochasticity of flux-freezing in not due intrin- 
sically to resistivity or other microscopic plasma mech- 
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Fig. 3. — Lagrangian trajectories that start in the ball of ra- 
dius p around space point x at time t move backward in time to 
to, explosively separating to a field-perpendicular distance Ax^ ~ 
{£\t - toP)^/^ which is independent of p for |t - to| > (p^/e)^/^. 



anisms that "break" field-lines but is, instead, a funda- 
mental consequence of turbulent Richardson diffusion. 



Higher-order terms in the generalized Ohm's law (47) 
that do not appear in the ideal equation (48) will leaoto 
melding and merging of field-lines at scales < pi. How- 
ever, the above argument strongly suggests that these 
details of the microscopic plasma processes do not af- 
fect the dynamics at scales larger than pi. In some cases 
this can be shown more analytically by defining a suit- 
able "motion" of field-lines consistent with the induc- 
tion equation. For example, the formulation in section 
4.1 based on addition of a Brownian motion to the La- 



grangian particle dynamics, eq.(31), can be carried over 
to certain instances of the generalized Ohm's law. See 
in particular Eyink (2009) for the Hall MHD equations. 
This approach is used in Appendix B to argue that nei- 
ther the Hall effect nor Ohmic resistivity will have any 
significant influence on the inertial-range turbulence dy- 
namics at large enough scales. The Hall term, for exam- 
ple, does not affect the dynamics at scales much greater 
than 5i = pi/\^Pi, the ion skin depth. Unfortunately, it 
is difficult to extend this type of argument to all cases 
of the generalized Ohm's law, because it is not known 
how to define a "motion" of field-lines consistent with 
the induction equation for the general case. 

On the other hand, there is a different argument which 
applies in general and leads to the same conclusion that 
flux-freezing must be intrinsically stochastic in turbu- 
lent plasmas. While the "motion" of magnetic field-lines 
is a conventional and somewhat arbitrary concept, the 
motion of plasma is perfectly well-defined within the va- 
lidity of an MHD description. Plasma fluid moves with 
the bulk velocity u. Thus, field-lines may be tracked by 
"tagging" the lines with plasma fluid element s and then 
following these as Lagrangian fluid particles (Newcomb 
1958 Axford 1984). In the case of a smooth, laminar so- 
lution of the ideal MHD equations, this is unambiguous 
because of Alfven's theorem: two plasma particles which 
start on a certain field-line must share a field-line for all 
times. One can then, by convention, consider this as the 
"same" field line as the initial one. This approach fails 





Fig. 4. — Advective acceleration of magnetic field-line separation. 
Top: A pair of ions starting on the same field-line. Middle: Electric 
field in the frame of the plasma changes the connections of the ions 
to field-lines separated by distance pQ. Bottom: The ions are swept 
apart by fluid motions to a distance ^ ^ po in a time t^ that is 
independent of po . The original field-line has effectively "split" into 
two separated by the macroscopic distance i. 



for a non-ideal Ohm's law. 



E 



1 



uxB = R, 



(51) 



where R represents all of the terms on the righthand side 
of (47) other than the motional term. Clearly, R is just 
the electric field in the rest frame of the plasma flow- 
ing with the bulk velocity u. Emf due to these non-ideal 
terms leads to time-dependent magnetic flux in the rest 
frame, corresponding to a slippage of field lines. This 
vitiates the usual method to assign an identity to in- 
dividual field-lines over time, because plasma elements 
shift their attachments to lines. Charged particles move 
along magnetic field-lines, but two plasma elements that 
start on one field-line will sit on distinct field-lines at 
later times. 

Now consider a turbulent plasma where the non-ideal 
term is numerically "small" but the plasma has a turbu- 
lent inertial-range in which the velocity field u is rough, 
with a power-law energy spectrum extending down to a 
smallest length-scale po ~ pi. The slight shifts in line- 
attachments are enormously amplified by explosive rela- 
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tive advection, as illustrated in Fig. 4. Consider a sin- 
gle magnetic field-line in this plasma and, along it, two 
plasma fluid particles at initial locations a and a^ Due 
to a combination of the non-ideal field R and advec- 
tion by sub-inertial-range eddies, the two plasma par- 
ticles will end up on distinct field-lines displaced a dis- 
tance |x(a',t) — x(a, t)| = po apart in a time tq which is 
generally microscopically small compared with the eddy- 
turnover time tL. Because of (44), the two plasma ele- 



ments will subsequently end up on field-lines displaced a 
macroscopic distance (et^)^/^ apart independent of po in 
a time t ^ Toj^ effect, the original field-line has split 
apart into twolines separated by macroscopic distance 
in a finite time t. Of course, this same argument works if 
we take a' = a + p for some small random displacement 
p along the initial field-line with \p\ < po . We have 
already argued in section 4.2 that the random ensemble 
of Lagrangian trajectories x(a + p, t) will coincide for 
t ^ Tq with the random ensemble x(a, t) obtained from 
the perturbation by Brownian motion in eq.(19). The 
initial field-line will be stochastically "frozen-in" to this 
ensemble of random flows. 

Turbulent advection accelerates the separation of field- 
lines so effectively that the microscopic plasma process of 
line-slippage, whatever its origin, is rapidly "forgotten" . 
The non-ideal terms in the generalized Ohm's law be- 
come insignificant in further separating lines which have 
exceeded the distance po apart. When the plasma flow 
velocity u is "rough" , it is the motional term in the gen- 
eralized Ohm's law which, after a very short time, domi- 
nates all of the other terms in the separation of field-lines. 



5. LARGE-SCALE MAGNETIC RECONNECTION 

With the theory developed above, we now turn to our 
main topic of magnetic reconnection. To be specific, we 
consider the steady-state reconnection of a pair of large- 
scale magnetic flux tubes which collide along a section 
of length Lx- We assume incompressible plasma flow and 
anti-parallel magnetic flux-tubes of equal field-strength 
\Bx\' The situation could be generalized to allow for com- 
pressible plas ma and asymmetric field-strengths ( [Cassak 



fc Shay^QQ7 ), flux-tubes intersecting at an angle or with 
a shared component of magnetic field (guide field) Bz 
twisted flux-tubes with magnetic 
All of these are 



1I2QQTI), 



etc. 



([Linton et al. 

helicity (Zweibel & R hoads | |1995) 
relevant for applications, but we keep the set-up simple 
in order to make the new ideas clear. The large-scale ge- 
ometry we consider is thus precisely the two-dimensional 
configuration in classical Sweet-Parker theory. However, 
we are working in three-space dimensions and thus small- 
scale turbulence, in particular, will be fully 3D. 

5.1. A Rederivation of LV99 by Stochastic 
Flux-Freezing 

The constraints on steady-state reconnection from the 
MHD balance equations are well-known. For a reconnec- 
tion layer of length and thickness A, the constraint 

-'^-'^ In the conventional view, for ideal plasmas with R = 0, the 
two elements would separate also to this distance but would remain 
on the same, highly-stretched field-line. This is not the case with 
a nonideality R 0, however small. 




Sweet -Parker model 





Turbulent model 



Fig. 5. — Upper plot: Sweet-Parker model of reconnection. The 
outflow is limited by a thin slot A, which is determined by Ohmic 
diffusivity. The other scale is an astrophysical scale ^ A. 
Lower plot: Reconnection of weakly stochastic magnetic fleld ac- 
cording to LV99. The model that accounts for the stochasticity 
of magnetic fleld lines. The outflow is limited by the wandering 
of magnetic fleld lines, which depends on fleld line stochasticity. 
From Lazarian et al. 2004. 



of energy conservation is 



1 

8^' 



(52) 



where Vrec is the in-flow plasma speed or reconnection 
velocity and Vout is the out-flow velocity from the sides 
of the reconnection layer. Here it is assumed that the 
energy in the out-flow jets is primarily kinetic, since the 
By component of the outgoing field-lines at the outer 
edges of the reconnection layer is still relatively small. 
The critical constraint of mass conservation is 



Vout^ 



Combining (52) and (53) yields 



Bx 



Vout 



(53) 



(54) 



SO that the out-flow velocity is at the Alfven speed of the 
in-coming magnetic fields. The momentum equation en- 
forces pressure balance along and across the reconnection 
layer and does not yield a new constraint in our setting. 
From (53) and (54) we see that 



(55) 



Our discussion to this point has been completely general 
and no particular mechanism of line-slippage has been as- 
sumed. To proceed further, a specific reconnection mech- 
anism must be invoked which permits one to estimate the 
thickness A of the layer. 

The classical Sweet-Parker model assumes resistive re- 
connection in a laminar MHD solution. There are then 
several ways to estimate the thickness A, which are all 
consistent. The most common aproach is to use the 
steady-state Faraday's law, which implies the constancy 
in the vertical ^/-direction of the reconnection electric 
field Ez. Outside the reconnection layer Ez = VrecBx/c, 
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while inside Ez = rjJz = tjcBx/AttA. Equating these 
gives 

A 



(56) 



where A 



(55) and 
tnat 



= T](? is the magnetic diffusivity. Combining 
([56]) gives the well-known Sweet-Parker results 

/\ = L^/^, Vrec = VA/^, (57) 

where S = vaLx/X is the Lundquist number. Another 
way to obtain this result is based on standard flux- 
freezing ideas for a laminar plasma flow (Kulsrud, 2005). 
The mean-square vertical distance that a magnetic field- 
line can diffuse by resistivity in time t is 



{y\t)) ~ Xt. 



(58) 



Each reconnected field-line is the result of resistive gluing 
of field-lines that diffuse vertically across the reconnec- 
tion layer. At the same time, the field lines are advected 
out of the sides of the reconnection layer of length at 
a velocity of order va- Thus, the time that the lines can 
spend in the resistive layer is the Alfven crossing time 
tA = Lx/va- Thus, field lines can only be merged that 
are separated by a distance 



(59) 



We thus recover ( |57[ ). This is an important dynamical 
consistency check on the Sweet-Parker model of laminar 
resistive reconnection. 

For large Lu ndq uist numbers, the Sweet-Parker recon- 
nection speed (57) is too small to account for observed 



reconnection rates in most astrophysical settings. How- 
ever, the Sweet-Parker laminar MHD solution is not ap- 
propriate for turbulent plasmas at large S. For sim- 
plicity, we shall consider reconnection in the presence 
of a steady-state background turbulence with an inte- 
gral length-scale Li and rms velocity ul- We shall con- 
sider sub- and trans- Alfvenic turbulence with ul ^ va ot 
Ma < 1, under the assumption that the large-scale mag- 
netic flux tubes are strong and well-defined in the pres- 
ence of turbulence. In some situations the reconnection 
may be initially laminar and then turbulence initiated by 
instabilities triggered by small fluctuations from t he en 



vironment. Initial fast reconnection by Hall effect dShay 
eTal] [19981 |Wang et al.]|2QQQl |Birn et al.|[2QQT| [Drake 
2001) or other collisionless processes may itself provide 
the energy release to generate small-scale MHD turbu- 
lence. We shall consider only the final stage of steady- 
state reconnection in the presence of turbulence with the 
given characteristic length and velocity scales. 

The key issue is to estimate A. We shall give two dif- 
ferent methods of determining the width of the turbulent 
reconnection layer. 

Our first estimate is based on stochastic flux- 
freezing/SFF (section 4) and 2-particle turbulent diffu- 
sivity (section 3.2). We have seen that GS95 theory im- 
plies that the turbulent separation of pairs of "virtual" 
field-lines backward in time, in the ^-direction perpen- 
dicular to the mean magnetic field, is given by eqs.(23) 
and (l28l), or 



{y\t)) 



U <^t <tc 



ULLiM^t, t > tc 



(60) 



for the strong and weak MHD turbulence regimes, re- 
spectively, with tc = Li/vA- Reconnected field-lines that 
emerge from the turbulent reconnection layer are con- 
stituted from field-lines obtained by following "virtual" 
field lines backward in time over the Alfven crossing pe- 
riod tA = Lx/vA' We thus obtain A = {y'^{tA)) as the 
width of the turbulent reconnection layer. There are two 
cases to consider, depending upon whether tA < tc or 
tA > tc Using £ ~ u\/vALi^ we rewrite (60) as 



{Ll/Li)M\{t/tAf iit<t, 
LiL^M\(t/tA), iit>tc. 



(61) 



Since tc ^ tA exactly when Li ^ Lx^ respectively, we get 
finally that 



A = L^M^min 




(62) 



This is precisely the result of Lazarian & Vishniac ( 1999 ). 
The present derivation avoids many complications m 
their original argument due to the fact that "virtual" 
field lines pass directly through each other. LV99 had 
to consider the limit to the reconnection rate set by the 
speed with which reconnected magnetic field elements 
can pass through each other on the way out of the recon- 
nection zone, employing a plausible but heuristic self- 
consistency argument. This problem is avoided com- 
pletely and rigorously by using SFF (which is mathe- 
matically equivalent to the resistive MHD equations). 

The thickness A of the turbulent reconnection layer 
was originally estimated in LV99 by a geometric argu- 
ment, based upon the spontaneous stochasticity of the 
magnetic field-lines themselves. The plasma in the re- 
connection zone can only be expelled along the lines of 
the (real, physical) magnetic field. As we have seen in 
section 3.1, the lateral wandering as one moves a distance 
s along a field line is given by eqs.([8|, ([9|, or 



isM%, 



\i s < Li 
if s > Li. 



(63) 



Since the plasma must move a distance s ^ Lx along the 
field-line in order to escape the reconnection layer, we can 
estimate the thickness of the layer by A ~ ^J((^p\Lx)) . 



Using (63) we immediately recover (62). The agreement 



of these two arguments demonstrates some basic dynam- 
ical consistency of the LV99 theory. 

We obtain from both arguments the LV99 estimate for 
the reconnection speed 



= vaM\ min 




(64) 



Note that v trans = vaM\ is the amplitude of velocity 
fluctuations at the transition between weak and strong 
MHD regimes. This can be a sizable fraction of the 
Alfven velocity. Furthermore, (64) is 'JasV^ in the 
strictest sense of being independent of the microscopic 
resistivity. The physical mechanism of this independence 
is spontaneous stochasticity, both of the Lagrangian par- 
ticle trajectories and of the magnetic field lines. In both 
cases, initial separations by a microscopic resistive length 
are amplified to macroscopic distances as one traverses 



16 



Eyink, Lazarian & Vishniac 



along the wandering curves by an amount (of time or 
of path-length, respectively) which is independent of the 
microscopic electrical resistivity. 

The above discussion has approximated the reconnec- 
tion layer as a zone of homogeneous MHD turbulence, 
with a constant Alfvenic Mach number Ma = ul/va set 
by the reconnecting x-component of the upstream mag- 
netic field, va = 5a:/y47rp. In reality, there will be a 
"magnetic shear layer" with a continuous change of the 
mean magnetic field Bx{y) across the layer. Likewise, 
there are nontrivial profiles Uy{x^ y) and Ux{x^ y) of both 
incoming and outgoing velocities, which are k nown to 



have nontrivial effects on 2-particle dispersion ( Shen 
Yeung 1997 Pope 2002). An interesting refinement of 



the original LV99 theory would be to consider the ef- 
fect of variation of the mean velocity and magnetic fields 
across the reconnection layer. One plausible effect of 
magnetic shear will be to increase Ma within the layer, 
permitting greater lateral wandering of field-lines and 
increasing the thickness A. Less easy to guess are the 
effects of the change in orientations of magnetic fields 
for flux-tubes reconnecting at an angle of less than 180° 
([Linton et al. 2001). The simulation results of Kowal] 
et al. " pb09| on tubes crossed at various angles give sup- 
port to the idea that the effects are minimal, but further 
theoretical and numerical study of this would be desir- 
able. Another effect which has been neglected in our 
derivation above is the turbulent energy dissipation that 
exists in the reconnection layer. This can be estimated 
as pe/S.Lx with e = u\/vALi. If this is added into the 



energy balance ( 52 ) one obtains a result for the outflow 
velocity 

f4 



1 - 2M\ 



(65) 



which can be substantially reduced from the Alfven 
speed. By the mass conservation constraint (53), there 
is a similar reduction in the reconnection speea. How- 
ever, these effects are ignorable under the condition 
M\ <^ Li/Lx or, equivalently, A <C L^. This is gener- 
ally well-satisfied. Other modifications and refinements 
of the LV99 model can be envisaged and are worth pur- 
suing (see the opening paragraph of this section and also 
section 6.4). 

It is important to emphasize that the original LV99 
argument made no essential use of averaging over turbu- 
lent ensembles. The stochastic line- wandering which is 
the essence of their argument holds in every realization 
of the flow, at each instant of time. The "spontaneous 
stochasticity" of field lines is not a statistical result in the 
usual sense of turbulence theory and does not arise from 
ensemble- averaging. The only use of ensembles in LV99 
is to get a measure of the "typical" wandering distance 
A of the field lines (in an rms sense). If one looks at 
different ensemble members of the turbulent flow, or at 
different single-time snapshots of the steady reconnection 
state, then A will fluctuate. Thus, the reconnection rate 
will also fluctuate a considerable amount over ensemble 
members or over time. E.g. see Figs. 12-14 in Kow al et] 
al. (2009). But it will be "fast" in each realization and at 
each mstant of time, because the mass outflow constraint 
is lifted by the large wandering of field lines. 

Exactly the same remarks apply to the rederivation of 
LV99 using stochastic flux- freezing (SFF) and Richard- 



son diffusion. SFF is also not "stochastic" in the usual 
sense of turbulence theory. The new derivation uses 
ensemble- aver aging just to get an estimate of the "typ- 
ical" wandering of stochastic Lagrangian trajectories in 
an rms sense. Numerical studies of Richardson's law 
in hydrodynamic turbulence show that it is a very robust 
property. At most points in a turbulent flow, a "cloud" of 
stochastic particles spreads with a power-law very close 
to t^. There is some intermittency effect, with certain 
rarer, high-intensity points showing a much faster dis- 
persion or low-intensity a slower dispersion, but it is not 
a huge effect. Only the higher-order moments of the rel- 
ative separation distance are strongly affected (Boffetta] 
fc Sokolov |[2QQ2l ). 

'I'he LV9y theory, therefore, does not involve "turbu- 
lent resistivity" or "turbulent magnetic diffusivity" as 
this is usually understood. This is ordinarily meant to 
be an enhanced diffusivity experienced by the ensemble- 
averaged magnetic field (B). However, it is only if one 
assumes some scale-separation between the mean and the 
fluctuations that the effect of the fluctuations can be 



legitimatel y described as an enhanced diffusivity (Mof- 
|fatt ||1983|). In realistic turbulent flows, with no scale- 
separation, this phenomenological description as an ef- 
fective diffusivity can be wildly inaccurate. LV99 makes 
no appeal to such concepts and, indeed, never consid- 
ers the ensemble-average field (B) at all. Neither does 
Nature. The Sun does not average over an ensemble of 
coronal loops in order to get fast reconnection! 

5.2. Small- Scale Reconnection and MHD Turbulence 

The focus of this paper has been on the reconnection of 
large-scale, oppositely-directed magnetic flux-tubes in a 
turbulent plasma, since it is this phenomenon which has 
the most direct interest for astrophysics. On the other 
hand, a similar process of reconnection of local magnetic 
field-lines must take place at small-scales in homogeneous 
MHD turbulence, with or without a background mag- 
netic field. It has been argued that such small-scale mag- 
netic reconnection is an essential feature of MHD turbu- 



lence |Matthaeus fc Lamkin 1986; Lazarian fc Vishniac 
1999| . A substantial literature has since arisen on the 
relation between MHD tur bulence and magnetic recon 
nection (e.g. se e Lazarian ( 2005 [ ~ ' ' 
Servidio et aL|(2010); Mininm 



Eyink fc Aluie j2006); 



Pouquet im<)l). We 



cannot review this entire subject here, but it is impor- 
tant to discuss briefly the implications of "spontaneous 
stochasticity" for the problem. 

It was observed already in LV99 that nearly parallel 
field-lines in adjacent turbulent eddies must frequently 
intersect and reconnect according to GS95 theory. The 
local geometry is similar to that in large-scale reconnec- 
tion "writ small" , with the local mean magnetic field act- 
ing as a guide-field and the field-lines with anti-parallel 
transverse components undergoing reconnection. The 
crossing-points and knots in the local magnetic field- 
lines, if unresolved, would create topological obstruc- 
tions to the plasma flow which would render hydrody- 
namic motion impossible. LV99 argued instead that 
line-intersections at perpendicular length-scale £± are re- 
solved by magnetic reconnection in the eddy-turnover 
time i±/5u{i)^ independent of the resistivity. Fast mag- 
netic reconnection is thus crucial for the very existence 
of MHD turbulence of the kind supposed in GS95. 
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The stochastic "line-wandering" implied by GS95 the- 
ory (section 3.1) suggests that a very complex line- 
topology must exist in turbulent flows with magnetic 
fields that are non-smooth on inertial-range scales. This 
complexity is made particularly evident by the "sponta- 
neous stochasticity" phenomenon (section 4.2), accord- 
ing to which there are infinitely-many field-lines passing 
though each point in the limit of zero resistivity! This 
result suggests that reconnection sites are distributed 
densely in space for MHD turbulence at very high con- 
ductivity. Turbulent plasma flows with rough magnetic 
and velocity fields have a very strange and intricate ge- 
ometry when compared with the smooth, laminar flows 
that have often been the focus of theoretical and em- 
pirical studies of reconnection. Another indication of 
this complexity is the fractal distribution of clustered 
magnetic nulls for a magnetic field with a turbulent, 
power-law energy spectrum (Albright 1999). This lat- 
ter fact is very intimately connected with the stochastic 
"line- wandering" phenomenon deduced in LV99, since 
the magnetic nulls form natural site s of separation of 
adjacent magnetic field-lines. See [Davila fc Vassilicos] 
[1 POOSI; 'Goto & Vassilicos | P004p for a discussion of 
the closely related hydrodynamic phenomenon associ- 
ated with Lagrangian particle trajectories and stagnation 
points of the turbulent velocity field. 

The above facts strongly suggest that magnetic re- 
connection must be a process that is active always and 
everywhere in homogeneous MHD turbulence. As one 
measure of the local reconnection rate one may take the 
estimates (23),(28) of the lateral diffusion of magnetic 
field-lines (Tor strong and weak turbulence, resp.) For 
example, in the strong turbulence regime, field lines dif- 



the coarse-grained velocity satisfies 
d 
dt 



A^'dx = c (j) s^'dx. (67) 



fuse in an eddy-turnover-time 



-^/^f/^ at scale 



i± by a mean-square perpendicular distance er/ = 
Thus, reconnection rates are exactly large enough to re- 
solve knots of field-lines at length-scale £± in the natural 
evolution time at that scale, as argued in LV99. Note 
that di±/dt Su{i) and thus the reconnection speed of 
pairs of field-lines separated by inertial-scale separations 
£± is equal to the relative velocity du{£) between the two 
lines. 

Another way to quantify small-scale reconnection in 
MHD turbulence is by violations of Alfven's Theore m at 
length-scale £. T his was the approach adopted by Eyink| 
& Aluie ( |2QQ6| ), who studied the MHD induction equa- 
tion m a coarse-grained description at the (isotropic) 
length-scale 4^ 

dtBi = Vx [u^xB^ + cs£ - AVxB^]. (66) 

Simple estimates show that the direct resistive contri- 
bution proportional to A is negligible at high magnetic 
Reynolds numbers whene ver the m agnetic energy stays 
finite in that limit. See Eyink | ( |2QQ 5). It is the turbulent 

subscale EMF si = (uxB)^ — u^xB^ which provides the 
necessary electric field for reconnection at the scale £. The 
magnetic flux through a material loop C£{t) advected by 

The field-perpendicular variations therefore include contribu- 
tions from eddies with down to £ but the field-parallel variations 
get contributions from somewhat smaller eddies with i± down to 



Eyink & Aluie (2006) established necessary conditions 
tor a non- vanishing line- volt age on the righthand side of 
the above equation, independent of £. These conditions 
are: either (i) a non-rectifiable (fractal) loop Ci{t)^ or (ii) 
a blow-up of u or B along the loop, or (iii) tangential dis- 
continuities of both u and B intersecting the loop Ci(t) 
over a finite length. One or more of these conditions is 
plausibly satisfied in turbulent flow, implying ubiquitous 
reconnection. An important remark is that MHD turbu- 
lence is dominated by sca le- local i nteractions under very 
general assumptions. See Eyink (2005); Aluie & Eyink 
(|2010j). In particular, the turbulent EMi^ is dominated 
by contributions from velocity and magnetic field modes 
around the scale £. It is the EMF from these modes which 
is responsible for reconnection of the lines of the coarse- 
grained magnetic fields at length-scale £ and microscopic 
plasma processes at smaller scales play no direct role. 

5.3. 2D MHD Turbulence and Reconnection 

While astrophysical reality is three-dimensional (3D), 
there is a scholarly and intellectual value to considering 
MHD turbulence in two-dimensional (2D) space. Com- 
parison of toy models with reality is often a good way to 
test and sharpen one's understanding. Numerical simu- 
lations can also be better resolved in 2D than in 3D, per- 
mitting study of higher kinetic and magnetic Reynolds- 
number flows. Furthermore, the conventional view is 
that MHD turbulence in 2D and in 3D are very similar, 
much more so than for hydrodynamics of neutral flu- 
ids. 3D MHD possesses the three ideal, quadratic invari- 
ants of total energy, cross-helicity and magnetic helicity 
and 2D MHD has a nearly identical set of ideal invari- 
ants, with magnetic helicity replaced by the mean-square 
magnetic potential. Absolute equilibrium spectra in 3D 
suggest that energy and cross-helicity cascade to high- 
wavenumb ers and m agnetic helicity to low-wavenumbers 
(Frisch et aL ||1975), while the same argument in 2D 
leads to an identical conclusion, with an inverse cascade 
of magnetic potential repla cing that of magnetic helicity 
( Fyfe fc Montgomery [1976 ). Numerical simulations have 
corroborated these predictions, with the most recent 2D 
MHD simulations giving energy spec tra close to k~y^ or 
in the forward cascade range dGomez et al. Il999 



Biskamp fc Schwarz |2001[|Ng et al 2003| ). On the other 
hand, from the point ot view ot G!Sy5 theory, 2D MHD 
turbulence appears as a very degenerate case. In partic- 
ular, shear-Alfven waves that play the dominant role in 
3D MHD turbulence according to GS95 are entirely lack- 
ing in 2D, where only pseudo-Alfven wave modes exist. 
This fact suggests that the character of MHD turbulence 
in 2D must be quite different than in 3D and certainly 
not well-described by GS95. 

Because of the numerical advantages, many turbulent 
reconnection studies have also bee n performed in 2D, be- 

f inning with the important work of |Matthaeus fc Lamkin| 
( 1985[ 1986). T he numerical evidence is rath er contro- 
versial, however. Matthaeus fc Lamkin ( 1985 ) reported 
the reconnection of magnetic islands at rates nearly in- 
dependent of resistivity in the early stages of their simu- 



18 



Eyink, Lazarian & Vishniac 



lations. Servidio et al. (2010) have more recently made a 
similar study of Ohmic electric fields at X points in homo- 
geneous, decaying 2D MHD turbulence. This is really a 
case of small-scale magnetic reconnection, as considered 
in our previous section 5.3, and not directly relevant to 
the issue of reconnection of large-scale flux tubes. Two 
other numerical studies have recently be en made of large- 
scale m agn etic reconnection in 2D, b y |Loureiro et al. 
(|2009D and |Kulpa-Dybel et al.| ( |2010| , wh ich reach dif 
lerent conclusions. Un the one hand, [Loureiro et al. 
12009) had a better resolution but used periodic bound- 
ary conditions, which strongly constrain the ability to do 
averaging of the reconnection rate and the attainment of 
the steady state for reconnection. They inferred from 
their data that the 2D turbulent reconnection rate may 
be in dependent of resistivity. On the other hand, Kulpa- 
ipybel et al.| (j2010) used smaller data cubes but longer 
averaging, which is enabled by their outflow boundary 
conditions. They concluded that the reconnection does 
depend on resistivity and therefore is slow. The raw data 
of the two groups are actually rather similar and higher 
Lundquist-number simulations are probably necessary to 
resolve the matter. 

What does theory say? In particular, can LV99 or some 
LV99-like theory apply in 2Dp^ At first thought, one 
might doubt that LV99 is relevant at all since its central 
element — stochastic line-wandering — seems ruled out in 
2D on the simple grounds that magnetic field-lines may 
not cross. However, closer consideration shows that this 
argument is fallacious. Magnetic field-lines can cross in 
2D, at very special points, the magnetic nulls. Further- 
more, in a rough magnetic field, the magnetic nulls and 
null- null line s will form a very complex, fractal web (Al- 
bright 1999) which may permit quite unconstrained wan- 



dering ot held-lines. Indeed, the phenomenon of sponta- 
neous stochasticity of field-lines ought to occur in 2D 
as well as in 3D, if the magnetic field is rough. There 
are rigorous examples of this type in 2D, where integral 
curves of a non-smooth vector field are non-unique at 
every po int in space (e.g. see section II. 5 and Fig. 2 in 
iHartman (2002).) From the Lagrangian point of view of 
Richardson diffusion and stochastic flux-freezing, there is 
also little difference between 2D and 3D. The Richardson 
t^-law has been observed, for example, in the inverse- 
energy cascade range of 2D hydrodynamic turbulence 



(Celani fc Vergassola 2001 Goto fc Vassilicos 2004) 



and its analogue could be expected in 2D MHD turbu 
lence. Indeed, if velocity fields are sufficiently rough in 
2D MHD turbulence, then spontaneous stochasticity of 
Lagrangian trajectories ought to hold in 2D just as in 
3D, and an LV99-like theory ought to apply. 

The above arguments are, however, quite delicate. As a 
counterexample, we may consider weak MHD turbulence 
in 2D. The most important difference from 3D is that 
only pseudo-Alfven wave modes exist in 2D. Because the 
weak turbulence cascade preserves /cy = 1/Li while in- 
creasing kj_, the total wavevector is nearly perpendicular 
and the pseudo-Alfven polarization vector nearly parallel 

■"^^ In a completely different class are recent approaches which 
attribute fast reconnection in 2D MHD to formation and merging 
of plasmoids without the presence of convention al turbulence in 
the current sheet (Samtaney et al. 2009 Cass ak et al. ||2009[ 
|Bhattacharjee et al. ^2009( |Huang &: Bhattacharje'e |2U1U| ). I Jiese 
will be discussed later m section <6.'6. 



to the magnetic field. Our estimate from section 3.2 for 
the 2-particle diffusivity of weak MHD turbulence is thus 
reduced in the field-perpendicular direction by a factor 
of {i±/i\\)'^. The result is that 



(68) 



Solving for the perpendicular particle separation from 
dl\/dt = D±{i) gives an exponential growth propor- 
tional to the initial separation 



£'i{t) ^i(0)exp((const.)6:t/v^). 



(69) 



i.e., no spontaneous stochasticity! There will be no 2D 
fast reconnection if such a weak turbulence cascade per- 
sists down to dissipation scales. Naively, one would ex- 
pect this to be true in 2D from a simple time-scale ar- 
gument. The interaction time of a pair of pseudo-Alfven 
waves with parallel length-scale i\\ and fluctuation ampli- 
tude Su£ is i\\/Su£^ which will in general be longer than 
the Alfven wave-period i\\/vA- This argument suggests 
that there is no strong MHD turbulence at all in 2D and 
thus no fast reconnection. 

The only way to avoid this conclusion that we can see is 
that the 2D turbulence develops very singular structures 
or "tangential discontinuities" with Su£ ~ va^ indepen- 
dent of £ over the inertia l range. Note that the rigorous 
result of Eyink & Aluie ( 2006 ) for the 2D case requires 
such discontinuities and/oFlblowup in the magnitudes of 
the velocity or magnetic fields. There seems also to be 
support for this scenario from numerical simulations of 
2D MHD turbulence, which show very intermittent cur- 
rent and vortex sheets in a background of weaker tur- 
bulence. It is possible that this is one of those instances 
where the weak-turbulence cascade is unsta ble to the for- 
mation of localized, singular s tructures (cf. Majda et al. 



n( [TM7t ; [Rumpf et al. 



^(2009)). The crucial issue seems 

to be whether these intermittent structures occur suffi- 
ciently densely in spacetime. We shall not attempt to 
resolve this issue here, but only conclude that the exis- 
tence of strong MHD turbulence and fast turbulent re- 
connection in 2D are quite open questions. 

5.4. Alternative Ideas on the Role of Turbulence in 
Magnetic Reconnection 

As we mentioned in the introduction, various ideas how 
turbulence can increase the reconnection rate were dis- 
cussed as far back as 40 years ago. However, these ideas 
fell short of solving the problem. For instance, some 
papers have concentrated on the effects that turbulence 
induces on the microphysical level. In particular, "Speiser] 
(1970) showed that in collisionless plasmas the electron 
collision time should be replaced wit h the ele ctron re ten- 
tion time in the current sheet. Also IJacobson fc MosesI 



(1984) proposed that the current diffusivity should be 
modilied to include the diffusion of electrons across the 
mean field due to small scale stochasticity. All these ef- 
fects can only marginally change the reconnection rates. 

The closest predecessor to Ly99 w as the important 
work of jMatthaeus fc Lamkin] ( [T985[ [l986) (ML85,86). 
Those authors studied 2D magnetic reconnection in the 
presence of external turbulence, both theoretically and 
numerically. They emphasized the very close analo- 
gies between the magnetic reconnection layer at high 
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Lundquist numbers and homogeneous MHD turbulence. 
They also pointed out various turbulence mechanisms 
that would enhance reconnection rates, including mul- 
tiple X-points as reconnection sites and motional EMF 
of magnetic bubbles advecting out of the reconnection 
zone. However, ML85,86 did not understand the impor- 
tance of "spontaneous stochasticity" of field lines and 
of Lagrangian trajectories and they did not arrive at an 
analytical prediction for the reconnection speed. An en- 
hancement of the reconnection rate was reported in their 
numerical study, but the setup precluded the calculation 
of a long term average reconnection rate. A more re- 



c ent st udy along the approach in Matthaeus fc Lamkin 
[1 ( |1985[ ) is one in | Wat son et al. | ( |2QQ7p , where the ef- 
Tects of small scale turbulence on 2D reconnection were 
studied and no significant effects of turbulence on re- 
connection were found. However, this new work stud- 
ies fiow-driven merging in the super- Alfvenic limit, with 
very strong fiows ramming magnetic field together. In 
this situation, turbulence-induced X-points that might 
provide additional reconnection sites are usually rapidly 
ejected from the sheet. As discussed in the previous sec- 
tion, it is still a wide open question whether fast recon- 
nection can occur in 2D MHD. 

Some other numerical simulation papers have claimed 
to see evidence for fast reconnection in resistive MHD 
simulations. An important study of tearing instabil- 
ity of current sheets in the presence of background 2D 
turbulence and the formation of large- scale, long-lived 
magne tic islands has been performed in |Politano et al.~| 
(1989). They present evidence for "fast energy dissipa- 
tion" in 2D MHD turbulence and show that this result 
does not change a s they change the resolut ion. A more 



recent work of of Mininni fc Pouquet | ( 2009) provides ev- 
idence for "fast dissipation" also in 3D MHD turbulence. 
This phenomenon is consistent with the idea of fast re- 
connection, but cannot be treated as a direct evidence of 
the process. Evidently, fast energy dissipation and fast 
magnetic reconnection are rather different physical pro- 
cesses. A se ries of p a pers b^^Ga lsgaard and Nordlund, 
in particular Galsgaard & Nordlund (1997b), might also 
be interpreted as providing indirect evidence for fast re- 
connection. The authors noted that in their simulations 
they could not produce highly twisted magnetic fields. 
One of the interpretations of this finding could be the 
relaxation of magnetic field via reconnection. In this 
case, this observations could be related to the numerical 



finding of [Lapenta fc Bettarini (2011) which shows that 
reconnecting magnetic configurations spontaneously get 
chaotic and dissipate, which in its turn may be related to 
the predictions in LV99. However, in view of many un- 
certainties of the numerical studies, we are not confident 
of this connection. The hi ghest resolution simulations 



of Galsgaard fc Nordlund (1997b) were only 136 and 



at modest Lundquist numbers that could not permit a 
turbulent inert ial-range. 

A number of papers have attempted to treat turbulent 
magnetic reconnection by the formal methods of Mean- 
Field Electrodyna mics (MFE), e.g. Strauss (1988); Kim 
fc Diamond (2001|). Before we address these works in any 



detail, let us brieiiy discuss the relationships of LV99 it- 
self to MFE. LV99 theory as originally formulated does 
not consider ensemble- or time- averaged magnetic fields, 
does not employ MFE and, in fact, makes only very lim- 



ited use of probability and averaging. The key observa- 
tion to derive LV99 predictions using the exact formula 
(35) for the magnetic field — the "stochastic Lundquist 
formula" — is that it does not reduce to the conventional 
Lundquist formula as A ^ in a turbulent plasma. 
The typical lateral spread of Lagrangian trajectories that 
contribute in that limit is estimated from GS95 turbu- 
lence theory to be £\ ~ e\t — to\^ for t < Li/vA and 
to be ^ ULLiM^\t - to\ for t > Li/vA- This is a 
turbulent mixing effect but not in the usual sense of 
an enhanced "eddy-diffusivity" or "/3-effect" experienced 
by a mean magnetic field. Within the latter concept 
(see Parker (1979j)), magnetic fields are presumed to be 
mixed up passively by turbulence and the rate of recon- 
nection thereby accelerated. This solution is, however, 
not tenable for any dynamical important magnetic field 
which would resist bending by turbulent motions at small 
scales. The concept in its original formulation is there- 
fore not applicable to any astrophysical fields, apart from 
unrealistically weak fields, which are of marginal astro- 
physical importance, anyhow. 

Of course, it is always possible to define in a purely for- 
mal manner an "effective diffusivity" which, when substi- 
tuted into the MFE equations, will reproduce the predic- 
tions of any theory for the reconnection velocity Vrec and 
width A of the reconnection layer. All that one needs to 
do is set 

Lx 9 



A 



turh 



(70) 



for the desired Vrec and for the corresponding layer width 
A = Lx{vrec/'^A) imposcd by mass balance. In the case 
of LV99, this leads to a "turbulent magnetic diffusivity" 



Ax 



turh 
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(71) 



where the subscript "_L" denotes that this diffusivity is in 
the direction transverse to the background field. There 
is not necessarily any real physics associated with this 
purely formal definition of a diffusivity. 

In the weak turbulence regime for Lx > Li, however, 
this "turbulent magnetic diffusivity" reduces to the tur- 
bulent thermal diffusivity derived by Lazarian (20061). 
This is not an accident. As we have discussed in "de- 



riving eq.(60), the nearly straight field-lines in the weak 



turbulence regime experience a Brownian motion in the 
plane perpendicular to their direction, with a diffusion 
constant X±turb = u^LiM^, as in ( |7l| ). This implies a 
term in the induction equation for the mean field B(x, t) 
of the form X±turb^±'^' There is therefore a consistent 
MFE approach to recover the predictions of the LV99 
theory for reconnection in the weak turbulence regime, 
by means of an appropriate anisotropic turbulent mag- 
netic diffusivity. 

This is not the true for the strong turbulence regime 
with Lx < Li. Recall from the discussion around eq.(60) 
that field lines in strong turbulence do not separate oii- 
fusively in the transverse direction, as i\{t) ~ X±turbt , 
but instead undergo a Richardson diffusion i^it) ~ et^. 
This superdiffusive line motion cannot be consistently de- 
scribed by a Laplacian term in the MFE equations. Note 
also that strict accuracy of the MFE description requires 
a slow variation of mean fields on the integral scale Li 
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of the turbulence, whereas the width of the reconnec- 
tion zone according to (62) always satisfies A < L^. One 



would therefore not expect MFE with a simple Laplacian 
diffusion term to consistently and accurately describe 
the reconnection layer in strong turbulence predicted by 
LV99. 

With this discussion as backdrop, let us review some 
of the attempts to d iscuss magnetic reconne ction using 



MFE, beginning with Kim & Diamond 
the first place, it must be emphasized t 



2001 ) [KDOl]. In 

lat tney consider 



a rather unconventional set-up with a very strong guide 
field {Bz) and with fluctuating components Brms also 
much larger than the reconnecting component (B//), i.e. 
{Bz) ^ Brms ^ (iBj:/ 1). Because of the first assumption, 
reconnecting lines of the mean field are almost exactly 
parallel, with only tiny components of opposing sign. Be- 
cause of the second assumption Ma ^ 1 and there is no 
recognizable reconnection geometry for actual lines of the 
unaveraged magnetic field: the anti-parallel horizontal 
components are just a small mean effect. These condi- 
tions were adopted for analytical convenience and do not 
generally reflect astrophysical reality (particularly the 
second). Then KDOl develop closure approximations for 
the turbulent magnetic diffusivity using the 3D RMHD 
equations. Assuming that the small-scale fluctuations 
are statistically stationary and appealing to boundary 
conditions (e.g. periodic) which permit one to neglect 
surface flux terms, they derive a strongly quenched tur- 
bulent diffusivity with Xturh = 0{X), the plasma diffusiv- 
ity. They therefore conclude that magnetic reconnection 
proceeds at the slow Sweet-Parker rate, in contradiction 
to LV99. However, we find their conclusions of little rel- 
evance to LV99. In the first place, they consider a quite 
different problem, with Ma ^ 1. In the second place, 
their main finding of a quenched turbulent magnetic dif- 
fusivity is due to the conservation of magnetic potential 
by the nonlinear terms of the RMHD equations. This 
quenching is unlikely to persist with the open boundary 
conditions required for stationary reconnection or to be 
a feature of turbulent reconnection governed by the full 
3D MHD equations. 

Finally, one should not mix up the concepts we dis- 
cuss here with the so-called "hyper-resistivity" ([Strauss 
1986! IBhattacharjee & Hameiri""1986( iHameiri <^ WsT 



tacharjee 1987; Diamond & Malkov 2003), which is an- 
other attempt to derive fast reconnection from turbu- 
lence within the context of mean- field resistive MHD. 
The form of the parallel electric field can be derived from 
magnetic helicity conservation. Integrating by parts one 
obtains a term which looks like an effective resistivity 
proportional to the magnetic helicity current. There are 
several assumptions implicit in this derivation. Funda- 
mental to the hyper-resistive approach is the assumption 
that the magnetic helicity of mean fields and of small- 
scale, statistically-stationary turbulent fields are sepa- 
rately conserved, up to tiny resistivity effects. However, 
this ignores magnetic helicity fluxes through open bound- 
aries, essential for stationary reconnection, that vitiate 
the conservation constraint. 

There is a general objection to all mean- field ap- 
proaches which, to "explain" fast reconnection, make 
appeal to some effective dissipation (resistivity, hyper- 
resistivity, etc.) experienced by the fields once averaged 
over ensembles or over small space- time scales. The dif- 



ficulty is that it is the lines of the full magnetic field that 
must be rapidly reconnected, not just the lines of the 
mean- field. The former implies the latter, but not con- 
versely. Nature does not average over ensembles of small- 
scale turbulence to get fast reconnection! No mean-field 
approach can claim to have explained the observed rapid 
pace of magnetic reconnection unless it is shown that the 
reconnection rates obtained in the theory are strictly in- 
dep endent of the length- and time-scales of the averaging. 
See Eyink fc Aluie | ([2006]) for further discussion. 



5.5. Turbulent Versus Collisionless Reconnection 

The most popular alternative to LV99 is currently not 
a competing turbulence theory of reconnection but in- 
stead non-MHD ap proaches based on the Hall effect in 
a two-fluid model (IShay et al.||1998l 119991 1 Wang et al." 
r^^f^. ■ -'^-^■ |Drake||Jb(jl| IM ' ' ' ' ' ' 



^2000! 'Birn et al. 2O011 1 

Cassak et al. 2010), or on fu 



alakit et al.||-iOa^ 

1 kinetic Vlasov dyna mics 



(Daughton et al. 20061 I^QOSl [2QTT| IChe et al. |201 1 ) 

have are;ued in section 4.3 and at length in Appenc 



We 

Appendix B 

that collisionless effects are, in most cases, irrelevant to 
determining reconnection speeds in the presence of tur- 
bulence. We do believe that there are some situations in 
certain astrophysical systems where kinetic effects will 
influence the rate of reconnection. However, formulating 
a criterion to determine if kinetic effects are important 
requires some care in the presence of turbulence. A pop- 
ular criterion for "collisionless reconnection" is the con- 
dition that the Sweet-Parker thickness S qp be smaller 



than the ion-skin depth 5i (Yamada et al. 2006 



Yamada 



Zweibel 



2009| ). This condition has been checked to 

correctly signal "collisionless reconnection" in laminar 
simulations and experiments. However, as we now dis- 
cuss, this condition in turbulent plasmas is not relevant 
to the validity of LV99. 

There is even some question of how to define the 
"Sweet-Parker thickness" 5sp in turbulent reconnection. 
In the original resistive MHD model of LV99, local, small- 
scale reconnection events were argued to have a Sweet- 
Parker Y-type structure and thickness 



L,M-'S-'^\ 



(72) 



with i:}^ the resistive cut-off length of GS95 turbulence 
theory. See LV99, eq.(14). However, this may not be 
the most relevant length-scale for determining the impor- 
tance of collisionless effects. Within GS95 phenomenol- 
ogy one can estimate the pointwise ratio of the Hall elec- 
tric field to the MHD motional field as 



J/en c6B{e^)/47rnee^ Si 1/2 
u Ul Li 



(73) 



where Sl = vaLi/X is the Lundquist number based on 
the forcing length-scale of the turbulence and Ma = 
ul/'^a is the Alfvenic Mach number. This suggests the 
definition of the "Sweet-Parker thickness" as 

5sp = LiMx'S-'/\ (74) 

SO that {J/en)/uL — Si/dsp- The length-scale Ssp is 
the analogue of the Taylor microscale in hydrodynamic 
turbulence, whereas the distance is the GS95 ana- 
logue of the Kolmogorov scale. If the magnetic diffu- 
sivity in the definition of the Lundquist number is as- 
sumed to be that based on the Spitzer resistivity, given 
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by A = 5lvth,e/^ei where 5e is the electron skin depth, 
Vth,e is the electron thermal velocity, and l^i is the elec- 
tron me an- free- path length for collisions with ions, then 

1/2 / . \ 2 



the plasma beta. Substituting into (73) and defining 
^ — '^th.i/'^* as the Mach number gives 



f m 

5sp ~ \me 



1/4 



M(3 



1/4 



Li 



1/2 



(75) 



which coinc ides precisely with the ratio defined by |Ya-| 
mada et al.| (2006 ) ,eq. (6) . 

With this aetinition of Ssp, it is easy to see from the 
data in our Table 1 that the ratio Si/Ssp — 10~^ for the 
warm ISM, 1 for post-CME current sheets and 10^ 
for solar wind impinging on the magnet osphere. How- 
ever, this has nothing to do with the validity of LV99 
for those systems. What it does say is that Ohmic re- 
sistivity may not be the mechanism for line-breaking in 
the latter two cases and the original resistive model of 
LV99 is not an adequate description at sufficiently small 
scales. Thus, the structure of local, small-scale reconnec- 
tion events should be strongly modified by Hall or other 
collisionless effects, with possibly an X-type structure, 
an ion layer thickness ~ Si, quadrupolar magnetic fields, 
etc. However, these local effects are irrelevant to deter- 
mining the global rate of large-scale reconnection (see 
Appendix B). 

What is the correct criterion to determine a break- 
down of LV99 and the importance of collisionless ef- 
fects? The LV99 model assumes that the thickness A 
of the reconnection layer is set by turbulent MHD dy- 
namics (line- wandering and Richardson diffusion). Thus, 
self-consistency requires that the length-scale A must lie 
within the range of scales where shear- Alfven modes are 
correctly described by incompressible MHD. This implies 
a criterion for "turbulent collisionless reconnection" 



(76) 

with A calculated from eq.(62) and pi the ion cyclotron 
radius (see section 2). Since~A oc L^, the large length- 
scale of the reconnecting ffux structures, this criterion is 
far from being satisfied in most astrophysical settings. 
For example, in the three cases of Table 1, one finds us- 
ing A = LM\ that Pi /A 10"^^ for the warm ISM, 
10~^ for post-CME sheets, and .1 for the magneto- 
sphere. Only in the latter case is the condition (76) close 



to being satisfied. Reconnection events in the magneto- 
sphere probably do typically involve collisionless effects 
in an essential way. It must be appreciated, however, 
that this is a highly nongeneric situation and turbulent 
ffuid effects will generally predominate in astrophysical 
reconnection. 

Even in magnetospheric reconnection the presence of 
background turbulence in the environment should not be 
ignored. Like nearly every cosmic plasma, the magneto- 
sphere commonly exhibits tur bulence, e ither statistically 
steady or in sporadic bursts (Zimbardo et al. 2010). En- 
ergy spectra are generally similar to those in the solar 
wind, with spectral exponents close to -5/3 and -7/3 at 
scales above and below the ion gyroradius/ion skin-depth 
(since /3 ^ 0.1 — 10, these are nea rly the sam e), respec- 
tively. See Table 1 in , Zimbardo et al] ( |2QTq1 ) . The orig- 



inal LV99 theory does not apply in the situation that 
Pi A > pe, but a version of LV99 based on EMHD 
might provide some insight. Of course, since pe is only 43 
times smaller than pi for a hydrogen plasma, any EMHD 
inertial range must be of limited extent. Kinetic effects 
enter at electron scales, not accurately described by a 
fluid model, and approaches based on solving the Vlasov 
equation (see section 6.3) become important. Some in- 
sights of LV99 should still carry over, such as the impor- 
tance of fleld-line wandering in enhancing reconnection 
rates. See section 6.3 for more discussion. 

5.6. Observational and Numerical Tests 

Let us discuss briefly the support for LV99 from obser- 
vations and simulations. 

The solar corona is a system proximate to Earth where 
LV99 theory should be applicable. Recent observations of 
a greatly thickened current sheet associated with coronal 
mass ejections give broad support to LV9 9 ( Ciaravella & 



Raymond||2008[ [Bemporad et al. ||2008| ). l^br example, 
the width A or the reconnection zone from the formula 



(62) was found by Ciaravella & Raymond |2008) to be 
just 10 times smaller than the observed thickness. This 
is quite good agreement, given that that formula has an 
unknown prefactor (which probably depends upon global 
geometry) and that some of the inputs (particularly the 
turbulence integral scale Li) were unknown from obser- 
vations and had to be crudely estimated. LV99 predicted 
also the phenomenon of triggering of magnetic reconnec- 
tion. Indeed, as the reconnection speed depends on the 
level of magnetic fleld stochasticity, the stochasticity in- 
duced by Alfven wave packages can enhance the recon- 
nection spe ed. This i s consistent with the recent obser- 
vations by |Sych et al.| ([2009) who reported a phenomeno- 
logical relationship between oscillations in a sunspot and 
flaring energy release above an active region above the 
sunspot. The authors proposed that the pulsations in 
the flaring energy release are triggered by wave packages 
arising from sunspots. The LV99 mechanism presents a 
very natural explanation for the phenomenon. 

While the Sun presents one of the best studied exam- 
ples of magnetic reconnection, interstellar medium has 
the advantage of testing LV99 in a collisional environ- 
ment (see |Yamada| ( |2007| )). Similarly, collisiona l recon- 
ne ction was observed in the solar photosphere (Park et 
al.||2009 ). For both types of environments effects of par- 
tial ionization may be important]^ The generalization 
of the LV99 model for a p artially ionized gas is presented 
in Lazarian et al.| ( 2004[ ). The direct numerical testing 
of the regimes ^oi reconnection described in that paper 
requires the use of a two-fluid code and has not been 
done so far. Given the importance of magnetic flelds in 
interstellar medium, such a testing is very appealing. 

The most convincing tests of LV99 would be of its 
predicted scaling laws for the width (62) and for the 
reconnection velocity (64) in terms of the parameters 



UL^VA^Li, and Lx- These scalings are most easily checked 
in numerical MHD studies, which achieve only moderate 
Reynol ds numbers but w hich permit controlled experi- 

\/vALi 



ments. Kowal et al. ( 2009 ) used the relation e 



For the interstellar media some phases are partially ionized 
(see [Draine &: Lazarian| ( |1998| for a list of idealized phases and 
their parameters). 
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to rewrite formula (64) as 



vaLx 



1/2 



mm{Lx,Li}. 



(77) 



and tested the s cahngs in s and L j by means of numeri- 
cal simulations. 'Kowal et al. (2009) studied in particular 
the weak-turbulence regime with > investigating 
the predicted power-law scalings Vrec e^^'^Li. The de- 
pendence on e^/^ was confirmed, but the predicted linear 
relation ex Li was not well verified and replaced with a 

weaker dependence closer to L^^^ for larger Li < L^. 
This behavior is possibly associated with a crossover to 
the strong-turbulence regime with Li > L^. Note that 
the formula ( 77 ) for strong turbulence predicts indepen- 
dence of the reconnection rate from L^, except through 
the dependence on Li of the mean energy dissipation e. 

The independence of reconnection rate on Li for strong 
turbulence is quite convenient for observational tests of 
LV99, since Li is much harder to reliably estimate in 
natural reconnection phenomena than are va^L^ and e. 
Indeed, the source of turbulence may not be so easily 
identified or localized in wavenumber. One ap proach is 



to obtain ul from Doppler line-broadening (e.g. Bempo 
rad et al. (2008)) and then to estimate Li = u\/vas. Of 



course, it spectral information is available, then one can 
obtain Li directly from the peak of the energy spectrum. 
The mean energy dissipation rate e is a source term for 
plasma heating, which can be estimated from observa- 
tions of electromagnetic radiation. To make such an es- 
timation reliably requires not only a radiation model but 
also some understanding of where the energy from the 
cascade is deposited at scales sma ller than the ion gyro- 
radius ( Schekochihin et al. ||2009| . It may be preferable 
to estiniate e from coarse-grained measurements of the 
energy-fiux rate 11^, which is equal, on average, to the dis- 
sipation rate. For example, the fiux of magnetic energy 
is given by 11^ = J^-s^ where = (c/47r)VxB£ is the 

coarse-grained current and = (uxB)^ — u^xB^ is the 
subscale EMF. Because of the scale- locality of the MHD 
energy cascade it possible to develop very quantitatively 
accurate approximations to the subscale EMF which de- 
pend only upon the coarse-grained gradients V u/ and 



VB^. For discussion of the approximation, see Eyink 



( |2006D . 

Additional testing of magnetic reconnection can be 
produced through indirect means. Indeed, the process 
of reconnection may be responsible for a wide variety of 
processes. For instance, the LV99 model of reconnec- 



tion w as invoked to explain gamma ray bursts ( Lazarian 
Zhang & Yan 201^ acceleration ot cosmic 
L^ 



|et al . 2002 

rays (de Gouveia Dal Pino & Tjazarian"2003^ 'Lazarian 
[2005 ; Lazarian Opher 2009 ; Lazarian & Desiati 2010), 
removal ot magnetic field during star formation (Lazar- 



ian||2005 


Santos-Lima et al. 


2010) 


quantifiec 


L and elaborated, it is p 



into the magnetic reconnection that is driving them. 

6. DISCUSSION 

Let us discuss here briefiy the main results of this paper 
as well as some of their ramifications. 

6.1. Flux- freezing in astrophysics 



Plasma conductivity is high for most astrophysical cir- 
cumstances. It has therefore generally been assumed that 
"fiux-freezing" is an excellent approximation for astro- 
physical plasmas. The principle of "frozen-in" field-lines 
provides a powerful heuristic which allows simple, back- 
of-the-envelope estimates in place of full solutions (ana- 
lytical or numerical) of the MHD equations. As such, the 
"fiux-freezing" principle has been applied to gain insight 
into diverse processes, such as star formation, stellar col- 
lapse, magnetic dynamo, solar wind-magnetospheric in- 
teractions, etc. The predictions of this simple principle 
often accord very well with observations. Flux-freezing is 
used to explain, for example, the magnetization of white 
dwarfs and neutron stars, the low angular momentum of 
stars, and the spiral structure of lines of force in the solar 
wind. 

However, for every success of the "fiux-freezing" prin- 
ciple, there is an equally stark failure. If the standard 
"fiux-frozenness" property held to good accuracy, then 
topology changes of magnetic field could not occur at 
the very fast rates observed in solar fiares and coro- 
nal mass ejections. During star- formation, the magnetic 
pressure of in-falling field-lines would be so great as to 
prevent gravitational collapse altogether. The tangled 
line-structure in small-scale dynamos would quench the 
exponential growth of magnetic field. The most com- 
mon attempts to explain such contrary observations in- 
voke additional physical effects that violate conventional 
"fiux-freezing", e.g. additional terms in the generalized 
Ohm's law besides collisional resistivity. Unfortunately, 
most of these microscopic plasma effects seem to be too 
small to have much effect in astrophysical MHD plasmas 
at very high kinetic and magnetic Reynolds numbers. 

Turbulence is a natural suspect to explain apparent 
breakdowns of "fiux-freezing" in astrophysical environ- 
ments. Most attempts to use turbulence invoke it as 
just another effective non-ideal term in the generalized 
Ohm's law. Indeed, any averaging of Ohm's law (e.g. 
over ensembles, space or time) produces an additional 
"turbulent EMF" that leads to a possible resistivity- 
independent violation of "fiux-freezing" for the mean 
magnetic field. This is a correct observation as far as 
it goes, but not a fundamental solution to the problem. 
"Flux-frozenness" in the conventional sense must be vi- 
olated for the full magnetic field and not just for some 
averaged field. 

Our solution is quite different. We have argued that 
"fiux-freezing" in the standard sense is violated in a 
turbulent inertial range at high kinetic and magnetic 
Reynolds numbers, but continues to hold in a novel 
stochastic sense. This stochastic fiux-freezing law in tur- 
bulent plasmas is a consequence of the remarkable tur- 
bulence phenomenon of "spontaneous stochasticity" of 
Lagrangian particle trajectories. As we have shown in 
detail, it is this stochastic form of "fiux-freezing" which 
underlies the LV99 theory of fast magnetic reconnec- 
tion. In any turbulent astrophysical system, "stochastic 
fiux-freezing" is the natural replacement for conventional 
fiux freezing. We believe that it will provide a powerful 
heuristic tool to explain many astrophysical phenomena. 
To use this tool with confidence requires a sound un- 
derstanding of MHD plasma turbulence, in particular 
Lagrangian particle statistics. This is a clear focus of 
further theoretical, experimental and numerical work. 
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6.2. Line-wandering and LV99 

Spontaneous stochasticity is a property not only of 
Lagrangian trajectories, but also of magnetic field-lines 
themselves. The identity of "the" field-line passing 
through a point becomes blurred in a high magnetic- 
Reynolds-number turbulent plasma. Certainly, at any 
moment in time there are a unique set of field lines that 
describe the field, through each point. However, there 
are an infinite number of field lines that can be drawn 
through a small volume, and if we follow for some parallel 
distance along these field lines, then they diverge explo- 
sively away from one another. In a turbulent plasma 
with a rough magnetic field, the distance that these lines 
diverge apart is independent of the dimensions of this ini- 
tial volume, when the parallel distance travelled greatly 
exceeds this dimension. As a limiting situation, for van- 
ishingly small initial volumes, there are infinitely many 
field-lines that emerge and meander randomly out of one 
point! At sufficiently large scales, that is, above the dis- 
sipation scale, this line- wandering is independent of the 
value of resistivity, and, within limits, to the details of 
the microscopic physics (see, for e xample. Appendix B). 

This was the original insight of Lazarian fc Vishniac 
(1999). The important implication is that plasma, flow- 



ing a finite distance along magnetic fields lines, will dif- 
fuse away from the mean magnetic field direction. This 
effect can accelerate reconnection speeds not only in the 
presence of inertial-range turbulence, but whenever field- 
lines exhibit some disordered, random character. If we 
use line-wandering in a turbulent reconnection zone to 
calculate the width of the outflow region around the cur- 
rent sheet, then we find that the classic Sweet-Parker rate 
becomes the fast reconnection rate given in LV99, inde- 
pendent of resistivity. As we have shown in the present 
paper, this same reconnection rate can also be obtained 
by an alternative approach in which an ensemble of field 
lines move stochastically in time. This agreement follows 
from a formal equivalence between the two approaches. 

The precise quantitative details of magnetic line- 
wandering depend upon the scaling laws of MHD tur- 
bulence. We have adopted here the predictions of phe- 
nomenological GS95 theory, and these should surely be 
verified and possibly refined in future numerical and ex- 
perimental studies. However, a parameter study in LV99 
shows that the conclusion of fast reconnection in a 3D 
turbulent flow holds for a range of spectral slopes and 
spectrum anisotropics that encompass the existing sug- 
gestions of modifying the GS95 model. Our work neglects 
also some features of large scale reconnection. By treat- 
ing the reconnection zone as a typical part of a turbulent 
plasma we are neglecting the effects of a strongly sheared 
large scale field and the fast outflow along magnetic fields 
that will accompany reconnection events. These features 
are natural ly included in numerical simulations (Kowal 
|et al.||2009|) and do not seem to result in any significant 
modification of the reconnection rate, but further work 
along this line is certainly desirable. 

6.3. Other recent work on reconnection 

Since our conclusions are consistent with the work of 
LV99 it is worth considering recent work on alternative 
approaches to calculating reconnection rates. Over the 
last decade, more traditional approaches to reconnection 



have changed considerably. At the time of its introduc- 
tion, the models competing with LV99 were modifications 
of the single X-po int collisionless r econnection scheme 
first introduced by 'Petschek (|1964|). Those models had 
point- wise localized reconnection regions which were sta- 
bilized via plasma effects so that the outflow opened up 
on larger scales (see Fig. 5). Such configurations would 
be difficult to realize in the presence of random forcing, 
which would be expected to collapse th e recon nection 
layer. Moreover, Ciaravella & Raymond ([2008") argued 
that observations of solar flares were inconsistent with 
single X-point reconnection. 

In response to these objections, more recent models of 
collisionless reconnection have acquired several features 
in common with the LV99 model. In particular, they 
have moved to consideration of volume-filling reconnec- 
tion, (although it is not clear how this volume filling is 
achi eved in the presen ce of a single reconnection layer 
(see [Drake et al. (2006|)). While much of the discussion 
still centers around magnetic islands produced by recon- 
nection, in three dimensions these islands are expected to 
evolve into contracting 3D loops or ropes due to tear ing- 
type in stab il ities in electron-current layers (Daug hton| 
|et al.[ ( [20Q8| |201ip ). This is broadly similar to wEat 
is depicted m h\g. 5, at least in the sense of introduc- 
ing stochasticity to the reconnection zone. The 3D PIC 
simulation studies reported in these works should help, 
in particular, to interpret magnetospheric reconnection 
phenomena. However, although the simulations are de- 
scribed as "turbulent", they do not exhibit the inertial- 
range power-law spectra observed in the magnetosphere 
and do not take into account either the pre-existing tur- 
bulence found in many of its regions (due to temperature 
anisotropy, velocity shear, Kelvin-Helmholtz instability, 
etc.) or inertial-range turbulence generated as a conse- 
quence of reconnection itself (Zimbardo et al. 201Q|. An 
EMHD analogue of LV99 an d^Kowal et al., ( ,2QQ9, ) studies 
may help to estimate such effects. 

Similar remarks apply to the recent 3D PIC study of 
|Che et al. [ ([2011|) , which observes micro-turbulence in 
the electron current layer during reconnection. The au- 
thors identify the source of this "turbulence" as a fila- 
ment ation instability driven by current gradients, very 
similar to a related instability in the EMHD model. The 
key aim of this work was to identify the term in the gen- 
eralized Ohm's law which supplies the reconnection elec- 
tric field to break the "frozen-in" condition. However, 
this study ignores the ambient inertial-range turbulence 
observed in the magnetosphere and other astrophysical 
plasmas, which may strongly modify laminar instabili- 
ties. Also, while there is interest in understanding the 
origin and character of reconnection electric fields (e.g. 
for particle acceleration) , we have argued at length in sec- 
tion 4.3 that the precise mechanism of "line-breaking" is 
irrelevant for the rate of reconnection in the presence of 
high- Reynolds- number inertial-range turbulence. 

Departure from the concept of laminar reconnection 
and the introduction of magnetic stochasticity is also ap- 
parent in a number of the recent papers appealing to the 
tearing mode instability to drive fast MHD rec onnection 



( [Loureiro et al. ||2QQ9[ [Bhattacharjee et al. |[2QQ9| ). [^ 



The idea of appealing to the tearing mode as a mea ns of 
enhancing the reconnection speed can be traced back to ,Strauss| 
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Several 2D numerical studies ( Samtaney et al. 2QQ9| 
Uassak et al. 2009; Bhattacha xiee et al |20091|Huan^ fc 
Bhattachariee 2010) have provided evidence that recon- 



,ttacharjee 2010) have provided 
iiection in the plasmoid-unstable region is independent of 
resistivity and a simple, heuris tic picture of a multi-level] 
plasmoid hierarchy has been proposed for this regime by! 
Uzdensk y'et al. | (|2010'rX recent very high-r esolution 
study of|i Ng Kagu nathan |(|2011|), however, h as found 
that plasmoids do not form even when the Eundquist 



five MHD mechanism of fast reconnection may require 



some level of evironmental noise, similar to LV99. In an 



case, since tearing modes exist even in a collisional fiuii 



this may open another chan nel of reconnection in such 
Muids. As we discuss below, this reconnection should not 



be "too fast" to account for the observational data. 



A f undamental consideration for all reconnection mod- 
els is that they must explain fast reconnection in colli- 



5y 

sional and collisionless plasmas. At the same time it| 



should be possible for reconnection to be delayed for sig- 
nificant amounts of time. Otherwise the accumulation 



of magnetic flux prior to a solar flare would be inexpli- 



cable. Fast reconnection as a feature of turbulence setsi 



number exceeds the putative critical value, except when 
insuffcient numerical resolution is used or when a small 



amount of noise is added externally. Thus, this alterna- 



ay 

3; 



a minimum r econnection spe ed which allows th e topol- 
ogy of the m agnetic field to evolve on dynamic al time 



scales. Thus the LV99 mode l can explain the accumu- 
lation of flux provided that th ere are n o generically fast 
reconnection proc esses at worlj^*^! An alternative expla-, 



nation based on c ollisionless reco nnection in a laminar] 
enviro nment has been suggest ed by|Cass ak et al. I(l2006l)l 
and [ Uzdensky (2007). Howe ver, it is not clear whether 
this explanat ion will account also f or the large observed 
thick ness of the macroscopic current s heets ( Ciar avella. 
^ Raymond 2008; Bemporad et al. 



if local turbulence is driven by release of energy from 



the magnetic field, it may result in a runaway turbulent 



reconnection process which may be relevant to some nu- 
merical simulations(Lapenta 2008; Bettarini & Lapenta 
^009). Alternatively, if tearing modes begin by driving 
relatively slow reconnection, but by some process desta- 



bilize the current sheet, then a similar runaway might re- 



sult. Turbulence looks like a natural candidate for such' 



process, but one should be open to a lternative sugges 
tions. For instance, we note that whil e X-point recon 
nection is clearly unstable for an isolated current sheet in. 



a collisional fiuid, interacting current sheets can produce 
bursts of fast X-point reconnection, separated by periods 
slow evolution ([Pang, Pen, fc Vishniac |(|2010[)). For- 



J ( [l988| >; |Waelbroeck| {1989}; [Shibata Tanuma | ( [200T]> . LV99 
showed tJiat tJie linear growtli ot tearing modes is msumcient to 
obtain fast reconnection. More recent work is based on the idea 
that the non-hnear growth of magnetic islands or plasmoids due to 
mergers provides large scale growth rates larger than the tearing 
mode linear growth rates on these scales. A situation where the 
non-linear growth is faster than the linear one is rather unu sual 
and requires further investigation. See Diamond et~al. | ( |198'4| . 

To be specific, take typical coronal loop paranieters n = 
10^ cm-3, T = 100 eV, B = 300 G, ul = lO'^ cm/s and U ^ 
La: = 10^ cm. Then va = 10^° cm/s and A = LiM\ = 10^ cm, 
comparable to the ion skin depth 6i = 700 cm. Thus, the initial 
turbulent reconnection rate will be quite small and may be slower 
than collisionless mechanisms. However, the energy released in 
the slow reconnection process can make the region more turbulent, 
accelerating the reconnection and resulting in a flare. 



mally, in this case we might not need to invoke turbulent 
feedbacl|^ but simply rely on the relatively slow pace of 
stochastic reconnection to evolve the magnetic field from 
one outburst to the next. 

In any case, in most astrophysical situations one has 
to deal with the pre-existing turbulence^ which is the in- 
evitable consequence of the high Reynolds number of as- 
trophysical fiuids and for which abundant empirical ev- 
idence exists. Such turbulence may modify or suppress 
instabilities, including the tearing mode instability. In 
this paper we have shown that it, by itself, induces fast 
reconnection on dynamical time scales. 

6.4. Outstanding Issues for Turbulent Reconnection 

Much remains to be done to fully clarify the effects 
of MHD plasma turbulence on astrophysical reconnec- 
tion. Without attempting to be exhaustive, let us men- 
tion some of the issues that we regard as most pressing. 

Of the possible extensions and refinements of LV99 
mentioned in section 5, we think one of the most im- 
portant is including the effects of magnetic and velocity 
shear. A related issue not addressed in the current pa- 
per is the detailed structure of the turbulent reconnection 
zone. We have invoked here the GS95 phenomenology of 
MHD turbulence, but intermittency effects not included 
in GS95 could be relevant. In principle, an independent 
theoretical derivation of the reconn ection speed can b e 
based on the line- voltage studied in Kowal et al. (2009). 
The arguments of LV99 suggest that the reconnection 
voltage should be contributed primarily by the motional 
EMF of already reconnected field-lines, but a detailed 
analysis is required. Addressing these issues will be the 
subject of future work. 

In its original formulation, LV99 is a model of non- 
relativistic reconnection. However, it is clear that the 
idea of enabling fast reconnection by extending the thick- 
ness of the outfiow region through magnetic field wander- 
ing should be applicable to relativist ic fiows. This exten- 
sion was i mplicitly used in some models of gamma ra y 



Zhang & Yan 



( |2QTT] )). 

ulent re- 



bursts (see Lazarian et al. (2002);^ 

Formal and detailed stumes of relativistic turbi 
connection are therefore important. 

Similarly, inertial-range turbulence with power-law 
spectra is an idealization adopted in LV99 to get ana- 
lytically tractable results. In many cases, however, wan- 
dering of field-lines and of Lagrangian trajectories is sig- 
nificantly affected by larger or smaller scales outside the 
inertial-range. Deviations from pure power-law spectra 
may occur due to multiple turbulent energy injection 
scales or to the effects of Ohmic dissipation. These ef- 
fects of inhomogeneity of turbulence driving are impor- 
tant subjects for further quantitative studies. 

Last, but not least, effects of turbulence dissipation 
via visc osity is another issue where more studies are re- 
quired. Lazarian et al. (2004) extended LV99 to the case 
of partially ionized gas where the dissipation of turbu- 
lence arises mostly through viscosity. The concept of 
Richardson diffusion and its relation to turbulent recon- 
nection in this case needs additional work. 

Energy release induces the feedback anyhow and the interac- 
tion of the resulting turbulence with the instabilities is an issue of 
further research. 
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7. SUMMARY 

The main points of our paper can be very briefly sum- 
marized as follows: 

1. In turbulent plasmas with rough velocity fields, the 
lines of force are not "frozen-in" in the standard de- 
terministic sense, but instead in a random sense as- 
sociated to the "spontaneous stochasticity" of La- 
grangian particle trajectories. 

2. Magnetic reconnection rates calculated based on 
"stochastic fiux-freezing" and phenomenological 
GS95 turbulence theory recover the predictions of 
the LV99 theory of turbulent reconnection. 

3. The LV99 predictions are independent of all non- 
ideal terms in the generalized Ohm's law. For 



example, they do not change if the Hall effect is 
included, whenever the ion skin depth is much 
smaller than the turbulent integral length or in- 
jection scale. 
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APPENDIX 

FORMAL TREATMENT OF TURBULENT LINE-SEPARATION 

Consider the problem of the growth in separation of a pair of magnetic field lines, starting at points displaced by 
vector as one moves a distance s in arclength along the field lines passing through the two points (cf. ,Jokipii | 
( 1973| )). The exact equation for the change in separation can be obtained from (45) to be 

±i{s) = h{^'{s))-h{m), (Al) 

where s is arclength, £{s) = ^\s) — ^(s), and b = B/|B| is the unit tangent vector along the magnetic field-line. By 
analogy with Richardson diffusion for Lagrangian particle trajectories, one can define a "2-line diffusivity" 

Dfj{£)= / ds{Sbi{£,0)Sbj{i,s)), (A2) 

^ — CO 

with units of length, where 6bi{i, s) = bi{^\s)) — bi{^{s)) with i = ^\0) — ^{0)^ so that 

^(€.(sK,(5)) = (AfW). (A3) 
The above equations are formally exact ( BatchelQ7]|1950[ [Kraichnan ||1966| ). Note that (A2) allows one to write 

Df^{£) 5k{i)Sb,{e)s,r^t{£) (A4) 

where Sint{^) is an integral correlation length of the increment in the tangent vector moving along the lines. (This is 
properly a tensor quantity, but here written as a scalar.) To determine the line wandering, one must have a model for 
this integral correlation length. 

The LV99 theory is obtained by modelling the integrand in ( |A2[ ) for field-perpendicular increments as 

{6b^{£,0)Sb^{£,s)) ^ ^!^Re(e^^/^ii-l^l/^(^)). (A5) 

Note that ~ SB{£)/Bo ~ Su{i)/vA' There are two effects in the s-dependent part. First, there is a correlation 

length A(^) of tangent- vector increments along the field line, which gives an exponential decay. Within LV99 theory, it 
is assumed that this length is the distance traveled by a random Alfvenic disturbance propagating along the field-line 
with velocity va in one energy cascade time at scale £±. That is, 

m = yAn = VA^-^. (A6) 

e 

The second effect arises from the periodic variation along the field-lines due to regular Alfven wave trains. It is assumed 
that, at perpendicular scale the Alfven wave trains that are dominant are those with wavelength ^y, the parallel 
length-scale corresponding to perpendicular scale i± . Integrating ( A5 ) in s gives the result 

Ai) - ^ iL - ^ JL (A7) 
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for A(^) > h. Substituting into (A4) gives 



Li 



A-> 



(A8) 



where notice that a factor of 5v?{t) has cancehed between numerator and denominator. From this last formula, various 
specific cases can be considered. 



In the strong GS95 turbulence regime, critical balance implies that \{€) ^ i\\ and thus, from (A7), Sint{^) ~ ^\\ - This 
gives the intuitive result that 

4/3 



M 



■4/3 
A •> 



(A9) 



where we have substituted from ([5|,(|6| for and Sui/va- Alternatively, one can just substitute for from ([5| into 
(A8). The equation 

4/3 

(AlO) 



is then easily solved to give eq.([8| in the text for Li > s ^ This result was already obtained for Ma = 1 by 
Skilling et al. ( 1974[) in a discussion of cosmic-ray diffusion, assuming a simple isotropic Kolmogorov hydrodynamic 
model ot magnetizea turbulence. In the weak turbulence regime one has instead h = Li Si constant, which substituted 



into (A8) gives 



ds 



LiM\ 



(All) 



and this is solved to give the diffusive result ^ in the text. 



NEGLECT OF THE HALL TERM IN TURBULENT RECONNECTION 



All of the nonideal terms in the generalized Ohm's law ( |47 ) 
-are expected 



-Ohmic resistivity. Hall term, pressure tensor, electron 



inertia, etc. — are expected to be insignificant in a long turbulent inertial range, at sufficiently large length-scales. 
Thus, they do not alter the predictions of the LV99 theory for the reconnection rates of large-scale flux structures. 
We shall here present more detailed analytical arguments for this thesis, in the specific example of the incompressible 
Hall MHD (HMHD) equations. As remarked in the text, the HMHD dynamics is valid in a literal sense for almost no 
astrophysical plasmas. Ho wever, it has been used in many reconnection studies as a simple fluid model that exemplifies 
th e Hall effect ( |Shay et al.|p98| p99| |Wang et al.l[2QQQ| |Birn et al.||2QQT| |Drake||2QQl| [Malakit et aL||2QQ9[ [Cassak et 
aL]|2 QlQ), so that it is a useful example to consider in this respect. 



he incompressible Hall MHD equations can be written as 



(StU + u-V)u = -Vp+jxb + i/Au, 
Stb = V X [(u - 5i]) X b] + AAb, 



(Bl) 
(B2) 



with V'U = V'b = and j = V xb. Here we have introduced the Alfven variable b = B/-y/47rp with units of velocity 
and 5i = c/up^i is the ion skin depth. The above system has two inviscid constants of motion that cascade to small 
scales, the total energy 

^ |2 



and the cross-helicity 



He 



E=\ yd3^[|u(x)|2 + |b(x)|2] 
= j (fxvi{-K) ■ [b(x) + i5iVxu(x) 



(B3) 



(B4) 



as well as the magnetic helicity which cascades to large scales. It is not our purpose here to give an exhaustive account 
of the turbulent cascade phenomenology of the HMHD model. Instead, we just wish to establish the simple fact that 
the Hall term modifies turbulence properties only at length-scales ^ ^ ^i, whereas length-scales ^ ^ will behave as 
in standard MHD turbulence without the Hall term. 



Our first argument is based on a spatial coarse-graining approach in an Eulerian formulation (e.g. see Eyink & Aluie 
[ ( 2QQ6[ )). The HMHD induction equation coarse-grained (low-pass filtered) at length-scale I becomes 



where the total contribution to the turbulent subscale EMF is 



and Up = u 



£l = (UeXb)^ - Ue^Xb^ 

6i} is the electron fluid velocity. It is easy to see that 



SiSb{£)/£ 



ul 



1/3 



(B5) 



(B6) 



(B7) 
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where the last expression uses GS95 theory to obtain the 1/3 powers of Ma and l/Li. However, a similar expression 
will hold in any theory of HMHD turbulence in which 5h{t) ^ VAM\{l/LiY with q>l and p > 0. Hence 5i]^ <C ui for 
<C ^ ^ I/i, so that the coarse-grained dynamics in HMHD turbulence at those scales is identical in form to standard 
MHD. The only difference is that the EMF (|B6|) has an additional contribution from the Hall term 



-J4(jxb)^ - j^xb^] = -5iV[ii{TMi)] 



T Mi 



(B8) 



where tm i = (bb)^ — b^b^ is the turbulent Maxwell stress. (Notice incidentally th at th e gradient term in (B8) does 



^MHU 



^Hall 



not contribute at all when substituted into the curl in (B5).) The total EMF in (B6) is thus 

including the standard MHD contribution s^^^ = (uxb)^ — u^xb^. It is easy to estimate the size of these various 
terms as 



^MHD 



5u{i)Sb{i), 



^Hall 



Si 



Sb\£) 



E 



■Ohm 



A 



5b{£) 



(B9) 



These estimates are rigorous upper bounds (cf. [Eyink (2005)) but they must be good order of magnitude estimates 
of the terms as well, in order to allow constant fluxes of energy and cross helicity to small scales. In that case we must 
interpret £ in the upper bounds as £±^ the scale perpendicular to the mean magnetic fielcp^ Instead at scales £ < 5i 
the terms proportional to Si dominate, so that 



n 



Hall 



Si 



Hall 



One thus gets constant fluxes with the scaling 

in agreement with the numerical results of [Shaikh fc Shukla] ( |2009 ) . As a matter of fact, the precise scaling at length- 
scales £ < Si will not affect our conclusions. As usual, the large-scale Ohmic electric field is negligible compared with 
the MHD turbulent EMF whenever £Su{£) > \ oi £ > £^ = {\^/ef^^ assuming GS9^^ However, one can also see 

that ef^^^ is negligible compared with e^^^ whenever £ ^ Si. Here we have assumed that Sh{£) Su{£)^ which is 
true in GS95 but also in any theory of MHD turbulence which predicts equipartition of kinetic and magnetic energies 
at small scales. These considerations show that the coarse-grained dynamics in Hall MHD turbulence reduce to those 
in standard MHD turbulence at length-scales £ ^ Si. This conclusion is in agreement with the numerical study of 
[Pmitruk fc Matthaeus (2006) . They compared results of Hall and standard MHD simulations at kinetic and magnetic 
Reynolds numbers of 1000, with Si for the Hall simulation chosen about 4 times greater than the resistive dissipation 
scale. They found virtually identical results in the two simulations at scales greater than Si for all of the various 
quantities, such as velocities, magnetic fields and electric fields. 

Now let us consider the Lagrangian description of HMHD turbulence. Stochastic flux-freezing holds also in resistive 
Hall MHD, with magnetic field-lines stochastically frozen-in to the electron fluid velocity Ue (Eyink ^^2009} . Thus, the 
lines can be regarded to "move" according to the stochastic equation 



u(x,t)-(5^j(x,t) + V2Ar7(t) 



(BIO) 



At separations £ ^ where the direct effects of the resistive term \/2Xri{t) can be neglected, the distance £±{t) = 
V that field-hnes advect apart in the field-perpendicular direction is obtained from 

d 



dt 



with 



Dm 



dt{5Ue±{£)-5Ue±{£,t)) 



(Bll) 



(B12) 



The electron velocity increment in (B12) gets a contribution Su{£) from the bulk plasma velocity and another contri- 
bution SiS}{i) from the Hall term. We would like to show that the Hall contribution is negligible whenever £ Si. 
This can be argued as follows: 



Note that energy flux to small scales is given in HMHD by 



- 

11^ - 

and cross-helicity flux by 



HV 



where Ct; = V X u is the vorticity and T£ = — tm £ is the total 
stress, both Maxwell stress and Reynolds stress tjh = (uu)^ — 



U£Ui. At scales £ > Si, the terms without 6i dominate and one gets 
constant energy and cross-helicity fluxes assuming GS95 scaling 
Su{i)r^6b{i) - {ei^y/^. 

For completeness we remark that the Ohmic fleld is also neg- 
ligible with respect to £^^^^ whenever 6i6b(i) > A. This deflnes 
another resistive scale i^"'^^ which is greater or smaller than Si 
according as £^ is greater or smaller than 6i. In the latter case, 
£Hall ± -g ^Yie true, physical resistive scale, not i:^ 
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First, we can estimate from energy balance alone without any other assumptions that e ~ ^jrms thus j ~ jrms = 
(6:/A)i/2. This allows us compare the typical size of the Hall velocity with the flow velocity as 

= (B13) 

u Ul ^ ^ ^ 

where we have used e = u\/vALi and we have defined = (A^/e)^/^, without the assumption that it plays any 
dynamical role as a "dissipation scale." Now let the Lundquist number Sl increase by increasing Li with e fixed, and 
thus also A,^^ and the ratio 

a = Si/£^ (B14) 



all fixed as well. It follows from (B13) that Sij <C u always when a < 1 but even when a > 1 if ^ Q^^. We 
therefore expect that the Hall contribution to the 2-particle diffusivity is negligible compared with the advective MHD 
contribution at high Lundquist numbers. 

Physically, the Hall term is very short-range correlated in both space and time, so that it acts on a pair of Lagrangian 
particles as two independent Brownian motions with an effective 1-particle diffusivity. This can be shown more formally, 
by observing that 

{Sh{£)-Sj±{£,t)) = 2[(ji(0)-ji(0,t)) - (j±{0)-h{£,m < 2(ji(0)-j^(0,t)), (B15) 

since the second term in the square brackets is a decaying term for increasing This upper bound is a considerable 
overestimate for smaller values of ^, but it suffices for our purpose. We can then estimate 





2 



dt{hiO)-hiO,t))^jtmsrx (B16) 

-CO 

where rx = (A/e)^/^ = l/jrms is the correlation time at the resistive dissipation scale. Multiplied by this is the Hall 
contribution to the 1-particle diffusivity. We therefore obtain a conservative upper bound on the Hall contribution to 
the 2-particle diffusivity, independent of ^, as 

Dl''"{^)<5ffr„.srx = a'X, (B17) 



where we have used (B13) and (B14). 

Finally, we compare with the advective contribution for i > S^. From the Eulerian argument presented earlier in the 
Appendix, we expect that GS95 scaling works on those length-scales. We can thus estimate for £j_ > Si that 

(5. 



4/3 

J^MHD^^^ ^ _ (^^4)1/3 _ ^4/3^ / ^ \ (g^g) 



where we have used A = {e£:^^y^^ and the definition of a. We conclude that 

Z)™'(£)/Z)f < (^] « 1 (B19) 



at least when ij_ ^ ^/o^Si. Thus, Richardson diffusion of field-lines in HMHD turbulence is exactly the same as for 
MHD turbulence at perpendicular separations i± much greater than the ion skin depth 5i. 

Now let us reconsider the problem of large-scale reconnection treated in section 5, including the Hall effect within 
the HMHD model. Define A = L^M| taking Li — L^. We shall make the essential assumption that (5^ ^ A. This 
assumption is almost always true in astrophysical reconnection problems. For example, using di = pij^fWi and the 
data in Table 1, we can see that 5i/ /S. ^ 10~^^ for the warm ISM, ~ 10~^ for post-CME current sheets and ~ 10~^ 
for solar wind impinging on the magnetosphere. Only in the latter case is the condition not extremely well satisfied. 
Let us then repeat the arguments in section 5 taking into account the Hall term. Consider first the derivation of LV99 
based on the stochastic wandering of magnetic field-lines as one traces points of increasing arc-length s along the lines. 
At small initial separations, the Hall term is important and the rms transverse distance = y^l^y^{s)) grows as 

i±{s) ^ {s'^ / 5iLi)M\ P] As soon as neighboring lines have separated to ^± ^ ^i, then their further separation is 
the same as in standam MHD turbulence without the Hall term. However, to reach the separation ~ 5i one must 

move along the field-lines a distance of only s ^ L^ = Li (^)^^^ <C Li when 6i <C A. In that case, the formula 



.A. 

,.2/ 



(63) in section 5.1 for {y {s)) still holds, if 5 L*. It therefore follows that the LV99 estimate for the width of the 



reconnection layer and the formula (64) for the reconnection velocity will hold, as long as 6i <C A, since most of the 



separation of lines will be achieved for s > L^. The same conclusion can be derived from the Lagrangian perspective. 

20 Note that 0-(r)jfe(O)>~A{56.^^^^^^^ 21 Obtained by solving M^/ds ^ dh{t)/vA 

tor large r wJienever t lie niagnetic-neld mcrenient ()b(r) lias scaling , , ^ \i/3 

exponent Kb < 1- {l/vA){ei^/Oi) . 
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< 5i the lines will separate laterally as i±{t) ~ (Jfet^)^/^ However, after a time 
{Li/va) much shorter than tc = Li/va the pairs of lines will reach the separation 



For small initial separations 

i±{t) ~ Si and GS95 scaling will take over. The formula (60) for {y'^{t)) that follows from GS95 will hold as soon as 
t > and the bulk of the separation of lines will occur after this time if Si <C A. Thus, the estimate (62) for the width 



of the reconnection layer holds just as before, with only a negligible correction due to finite Si 

We thus reach the important conclusion that the Hall term is irrelevant to determining the rate of global reconnection 
in the presence of high- Reynolds-number plasma turbulence. At most a modest enhancement of the reconne ction rate 
is obtained , but vanishingly small when Si <C A. We should mention that a similar conclusion was reached by [Smith et 



al. 



(2004) on the basis of numerical simulations. However, that study can be fairly criticized because the simulations 



were performed in 2D and at such low Reynolds numbers that most of the measured reconnection rates (in a transient 
regime) were smaller than the steady-state Sweet-Parker rates. Our theoretical arguments, on the other hand, are valid 
in 3D and at arbitrarily high Reynolds numbers. This is not to say that the Hall term will have no effect whatsoever, 
but these will be confined to very small scales, e.g. modifying the structure of local reconnection events. However, 
these local effects are irrelevant to determining the global rate of large-scale reconnection. 

Similar conclusions can be reached not only for the 2-fiuid Hall term, but also for other microscopic plasma mecha- 
nisms proposed to enhance reconnection speeds. For example, the MHD turbulent cascade may be terminated not by 
Spitzer resistivity but instead by some form of anomalous resistivity, e.g. due to ion acoustic instability. The latter 



produces a large e ffecti ve magn etic diffusivity of order {mi/rUeY^'^SluJc^e at scales smaller than Si/y/P (e.g. see Galeev 



& Sagdeev (1984) and Kulsrud ^ (2005 J , section 14.7.) But larger scales will behave as stand ard MHD turbulenc e and 



the LV99 predictions will again hold as long as (5^ <C A. This was verified in the simulations of [Kowal et al. (2009), who 
found no effect of anomalous resistivity on reconnection speed. Turbulence is such a powerful accelerator of magnetic 
line-separation that it completely dominates microscopic plasma processes of line diffusion, as soon as the magnetic 
and kinetic Reynolds numbers are sufficiently high. 



22 Obtained by solving M^/dt = Si6b{£)/£^ = {e6^/£^)^/^ 
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